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Abstract 

Atoms with a large scattering length have universal low-energy properties that do 
not depend on the details of their structure or their interactions at short distances. 
In the 2-atom sector, the universal properties are familiar and depend only on the 
scattering length. In the 3-atom sector for identical bosons, the universal properties 
include the existence of a sequence of shallow triatomic molecules called Efimov 
trimers and log-periodic dependence of scattering observables on the energy and 
the scattering length. In this review, we summarize the universal results that are 
currently known. We also summarize the experimental information that is currently 
available with an emphasis on 3-atom loss processes. 
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1 Introduction 



The scattering of two particles with sufficiently low kinetic energy is deter- 
mined by their S-wave scattering length, which is commonly denoted by a. The 
energies of the particles are sufficiently low if their de Broglie wavelengths are 
large compared to the range of the interaction. The scattering length a is also 
the most important interaction variable for 3-body systems and for higher 
A'"-body systems if all their constituents have sufficiently low energy. 

Generically, the scattering length a is comparable in magnitude to the range 
Vq of the interaction: \a\ ~ Tq. In exceptional cases, the scattering length can 
be much larger in magnitude than the range: \a\ 3> Tq. Such a large scattering 
length requires the fine-tuning of a parameter characterizing the interactions 
to the neighborhood of a critical value at which a diverges to ±cxd. If the 
scattering length is large, the atoms exhibit properties that depend on a but 
are insensitive to the range and other details of the short-range interaction. 
We will refer to these properties as universal, because they apply equally well 
to any nonrelativistic particle with short range interactions that produce a 
large scattering length. 

In the 2-body sector, the universal properties are simple but nontrivial. For 
example, in the case of equal-mass particles with a > 0, there is a 2-body 
bound state near the scattering threshold with binding energy Er, — h'^/ (ma^). 
The corrections to this formula are parametrically small: they are suppressed 
by powers of ro/a. 

In the 3-body sector, the universal properties are much more intricate. The 
first strong evidence for universal behavior in the 3-body system was the 
discovery of the Efimov effect by Vitaly Efimov in 1969 [1]. In the resonant 
limit a — > ±cxd, there is a 2-body bound state exactly at the 2-body scattering 
threshold E — Q. Efimov showed that in this limit there can also be infinitely 
many, arbitrarily-shallow 3-body bound states whose binding energies 
have an accumulation point at = 0. They are called Efimov states. As the 
threshold is approached, the ratio of the binding energies of successive Efimov 
states approaches a universal constant. In the case of identical bosons, the 
asymptotic ratio is 

E^^^^'^/E^^ — ^1/515.03 , as n ^ +00 with a = ±oo . (1) 

The universal ratio in Eq. (1) is independent of the mass or structure of the 
identical particles and independent of the form of their short-range interac- 
tions. The Efimov effect can also occur in other 3-body systems if at least 
two of the three pairs have a large S-wave scattering length. The numerical 
value of the asymptotic ratio may differ from the value in Eq. (1). Efimov's 
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discovery was at first greeted with some skepticism. However the theoretical 
evidence for the Efimov effect quickly became conclusive. 

The Efimov effect proved to be just the first nugget from a gold mine of 
universal aspects of the 3-body problem. This system has universal properties 
not only in the resonant limit a — ±00, but whenever the scattering length is 
large compared to the range tq. In two brilliant papers in 1971 and 1979 [2,3], 
Efimov derived a number of imiversal results on low-energy 3-body observables 
for three identical bosons. The dependence of these results on the scattering 
length or the energy is characterized by scaling behavior modulo coefficients 
that are log-periodic functions. This behavior is characteristic of a system with 
a discrete scaling symmetry. We will refer to universal aspects associated with 
a discrete scaling symmetry as Efimov physics. 

Although the existence of Efimov states quickly became well-established the- 
oretically, their experimental confirmation has proved to be more challenging. 
Although more than 36 years have elapsed since Efimov's discovery, there 
has still not been any convincing direct observation of an Efimov state. One 
promising system for observing Efimov states is ^He atoms, which have a 
scattering length that is more than a factor of 10 larger than the range of 
the interaction. Calculations using accurate potential models indicate that 
the system of three ^He atoms has two 3-body bound states or trimers. The 
ground-state trimer is interpreted by some (including the authors) as an Efi- 
mov state, and it has been observed in experiments involving the scattering 
of cold jets of ^He atoms from a diffraction grating [4]. The excited trimer is 
universally believed to be an Efimov state, but it has not yet been observed. 

The rapid development of the field of cold atom physics has opened up new 
opportunities for the experimental study of Efimov physics. This is made pos- 
sible by two separate technological developments. One is the technology for 
cooling atoms to the extremely low temperatures where Efimov physics plays 
a dramatic role. The other is the technology for controlling the interactions 
between atoms. By tuning the magnetic field to a Feshbach resonance, the 
scattering lengths of the atoms can be made arbitrarily large. Both of these 
technological developments were crucial in a recent experiment that provided 
the first indirect evidence for the existence of an Efimov state [5]. The sig- 
nature of the Efimov state was a resonant enhancement of the loss rate from 
3-body recombination in an ultracold gas of ^^^Cs atoms. 

This experiment is just the beginning of the study of Efimov physics in ul- 
tracold atoms. We have recently written a thorough review of universality in 
few-body physics with large scattering length [6]. In this shorter review, we 
summarize the results derived in Ref. [6], include a few more recent develop- 
ments, and focus more directly on applications in atomic physics. 
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In Section 2, we introduce various scattering concepts that play an important 
role in Efimov physics. Most of this review is focused on identical bosons, 
because this is the simplest case in which Efimov states arise and because it 
is the system for which Efimov physics has been most thoroughly explored. In 
Section 3 and 4, we summarize the universal features of identical bosons with 
a large scattering length in the 2-body and 3-body sectors, respectively. In 
Section 5, we describe the effects of deep diatomic molecules on the universal 
results. In Section 6, we describe applications to "^He atoms and to alkali atoms 
near a Feshbach resonance. In Section 7, we discuss the conditions under which 
Efimov physics arises in systems other than identical bosons. In Section 8, we 
consider the corrections to the universal results from the nonzero effective 
range and from microscopic models of the atoms as well as the extension of 
universality to 4-body systems. We conclude with the outlook for the study 
of Efimov physics in ultracold atoms. 



2 Scattering Concepts 

2.1 Scattering length 

The elastic scattering of two atoms of mass m and total kinetic energy E — 
h^k"^ /m can be described by a stationary wave function '0(r) that depends on 
the separation vector r of the two atoms. Its asymptotic behavior as r — > oo 
is the sum of a plane wave and an outgoing spherical wave: 

^{r) ^ e^^^ + h{e)— , (2) 

r 

where fk{0) is the scattering amplitude, which depends on the scattering an- 
gle 6 and the wave number k. The differential cross section da/dfl for the 
scattering of identical bosons can be expressed in the form 

^-\fk{o) + Mn-e)f . (3) 

If the two atoms are identical fermions, the + should be replaced by — . If they 
are distinguishable atoms, the term +/fc(7r — 6*) should be omitted. The elastic 
cross section a{E) is obtained by integrating over only | the Att solid angle if 
the two atoms are identical bosons or identical fermions and over the entire 
An solid angle if they are distinguishable. 

The partial-wave expansion resolves the scattering amplitude fk{0) into con- 
tributions from definite angular momentum quantum number L by expanding 



4 



it in terms of Legendre polynomials of cos^: 



oo 



2L + 1 



fk{e) = E 



k cot SL{k) — ik 



Pl(cos^), 



(4) 



L=0 



If there are no inelastic 2-body channels, the phase shifts dL{k) are real-valued. 
If there are inelastic channels, the phase shifts can be complex-valued with 
positive imaginary parts. 

If the atoms interact through a short-range 2-body potential, then the phase 
shift approaches zero like in the low-energy limit k 0. Thus 

S-wave (L = 0) scattering dominates in the low-energy limit unless the atoms 
are identical fermions, in which case P-wave (L = 1) scattering dominates. 
At sufficiently low energies, the S-wave phase shift So{k) can be expanded in 
powers of k'^ [7] . The expansion is called the effective-range expansion and is 
conventionally expressed in the form 



where a is the scattering length and the S-wave effective range. 
2.2 Natural low-energy length scale 

At sufficiently low energies, atoms behave like point particles with short-range 
interactions. The length scale that governs the quantum behavior of the center- 
of-mass coordinate of an atom is the de Broglie wavelength A = 27rh/p, where 
p is the momentum of the atom. If the relative momentum p of two atoms 
is sufficiently small, their de Broglie wavelengths are larger than the spacial 
extent of the atoms and they are unable to resolve each other's internal struc- 
ture. Their interactions will therefore be indistinguishable from those of point 
particles. If the atoms interact through a short-range potential with range 
To and if the relative momentum of the two atoms satisfies p <^ h/ro, then 
their dc Broglie wavelengths prevent them from resolving the structure of the 
potential. 

For real atoms, the potential is not quite short-range. The interatomic po- 
tential V{r) between two neutral atoms in their ground states consists of a 
short range potential and a long-range tail provided by the van der Waals 
interaction: 

V{r) — ^ as r — > oo . (6) 



k cot do{k) — —l/a-\- ^rgk'^ + . . . , 



(5) 
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The van der Waals tail of the potential does not prevent the scattering am- 
plitude from being expanded in powers of the relative momentum fc, but at 
^th Q];(^gj; i]^Q dependence on k becomes nonpolynomial. The coefficient Cq 
determines a length scale called the van der Waals length: 



At sufficiently low energy, the interactions between atoms are dominated by 
the van der Waals interaction. Thus the van der Waals length in Eq. (7) is the 
natural low-energy length scale for atoms. 

The natural low-energy length scale sets the natural scale for the coefficients 
in the low-energy expansion of the scattering amplitude fk{d)- It is sometimes 
referred to as the characteristic radius of interaction and often denoted by 
tq. If the magnitude \a\ of the scattering length is comparable to i, we say 
that a has a natural size. If |a| 3> we call the scattering length unnaturally 
large, or just large to be concise. The natural low-energy length scale also sets 
the natural scale for the effective range defined by Eq. (5). Even if a is 
large, we expect to have a natural magnitude of order i. For a and to 
both be unnaturally large would require the simultaneous fine-tuning of two 
parameters in the potential. 

The natural low-energy length scale £ also sets the natural scale for the binding 
energies of the 2-body bound states closest to threshold. The binding energy 
of the shallowest bound state is expected to be of order h^/mi'^ or larger. It 
can be orders of magnitude smaller only if there is a large positive scattering 
length a > 

2.3 Atoms with large scattering length 

The scattering length a can be orders of magnitude larger than the natural low- 
energy scale only if some parameter is tuned to near a critical value at which a 
diverges. This fine-tuning can be due to fortuitous values of the fundamental 
constants of nature, in which case we call it accidental fine-tuning, or it can be 
due to the adjustment of parameters that are under experimental control, in 
which case we call it experimental fine-tuning. We will give examples of atoms 
with both kinds of fine-tunings. 

The simplest example of an atom with a large positive scattering length is 
the helium atom ^Hc. The van der Waals length defined by Eq. (7) is i^dw ~ 
10.2 floi where Oq = 0.529177 A is the Bohr radius. The scattering length can be 
calculated precisely using potential models for helium atoms. For example, the 
LM2M2 [8] and TTY [9] potentials have a large scattering length a = 189 Oq 




(7) 
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but a natural effective range = 14 oq. Tliey predict tliat ^He atoms liave a 
single 2-body bound state or dimer, which is very weakly bound. The binding 
energy of the dimer is E2 = 1.31 mK, which is much smaller than the natural 
low-energy scale h'^/mi'l^^^ ^ 420 m The scattering length of ^He atoms 
is large because of an accidental fine-tuning. The mass of the "^He nucleus, the 
electron mass, and the fine structure constant a of QED have fortuitous values 
that make the potential between two ^He atoms just deep enough to have a 
bound state very close to threshold, and therefore a large scattering length. If 
one of the ^He atoms is replaced by a '^He atom, which decreases the reduced 
mass by 14% without changing the interaction potential, the scattering length 
has the more natural value — 33ao. 

The simplest example of an atom with a large negative scattering length 
is the polarized tritium atom [10]. The van der Waals length for is 
^vdw = 13.7 Oq. The scattering length for polarized '^H atoms is the spin- 
triplet scattering length at = — 82.1ao [10], which is much larger than £vdw- 
Polarized tritium atoms have no 2-body bound states, but they have a single 
3-body bound state with a shallow binding energy of about 4.59 mK [10]. 

Other examples of atoms with large scattering lengths due to accidental fine 
tuning can be found among the heavier alkali atoms. The spin-triplet scatter- 
ing lengths at for ^Li and for ^'^^Cs and the spin-singlet scattering length for 
^^Rb are all more than an order of magnitude larger than the corresponding 
van der Waals scales ^vdw- The fine-tuning is illustrated by the facts that ^Li, 
whose mass is 17% larger than that of ^Li, has a natural value for at and that 
*^Rb, whose mass is 2.3% larger than that of '^^Rb, has a natural value for a^. 

The mechanism for generating a large scattering length that involves tuning 
the depth or range of the potential is called a shape resonance. With this 
mechanism, only the open channel defined by the scattering particles plays an 
important role. Another mechanism for generating a large scattering length 
is a Feshbach resonance [11]. This requires a second closed channel in which 
scattering states are energetically forbidden that is weakly coupled to the 
open channel. A large scattering length for particles in the open channel can 
be generated by tuning the depth of the potential for the closed channel to 
bring one of its bound states close to the threshold for the open channel. The 
resulting enhancement of the scattering of particles in the open channel is a 
Feshbach resonance. 

Feshbach resonances in alkali atoms can be created by tuning the magnetic 
field [12,13]. In this case, the open channel consists of a pair of atoms in a 
specific hyperfine spin state \ f,mf). The closed channel consists of a pair of 
atoms in different hyperfine states with a higher scattering threshold. The 



The conversion formula for this energy unit is 1 mK = 8.61734 x 10 eV. 
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weak coupling between the channels is provided by the hyperfine interaction. 
Since different hyperfine states have different magnetic moments, a magnetic 
field can be used to vary the energy gap between the scattering thresholds and 
bring a bound state in the closed channel into resonance with the threshold 
of the open channel. The resulting enhancement of the scattering of particles 
in the open channel is a Feshbach resonance. The scattering lengths generally 
vary slowly with the magnetic field B. If there is a Feshbach resonance at 
Bj-es, the scattering length varies dramatically with the magnetic field in the 
vicinity of Bj-^g- If the Feshbach resonance is narrow, the scattering length near 
the resonance has the approximate form 



If B is far below or far above i?res, the scattering length a{B) has the off- 
resonant value abg. The parameter Ares, which controls the width of the reso- 
nance, is defined so that a{B) vanishes when B — B^^s + Ares- As B increases 
through -Brcs, a{B) increases or decreases to ±cxd, jumps discontinuously to 
=Foo, and then returns to its off-resonant value. The magnetic field provides 
an experimental fine-tuning parameter that can be used to make \a\ arbitrar- 
ily large. The existence of Feshbach resonances in atomic physics was first 
predicted in Ref. [13] for the specific case of Cs atoms. The use of a Feshbach 
resonance to produce a large scattering length in alkali atoms was first demon- 
strated by the MIT group using experiments with Bose-Einstein condensates 
of ^^Na atoms [14,15] and by the Texas and JILA groups using experiments 
with ultracold gases of ^^Rb atoms [16,17]. 

For ^"^Na atoms, the spin-singlet and spin-triplet scattering lengths are both 
smaller than the natural scale ^vdw, so all the hyperfine spin states \f,mf) have 
natural scattering lengths a,t B — 0. However, there are Feshbach resonances 
at which the scattering lengths diverge. For example, the ]1,+1) hyperfine 
state has Feshbach resonances near 853 G and 907 G [14] and the |1,— 1) 
hyperfine state has a Feshbach resonance near 1195 G [15]. The Feshbach 
resonances near 853 G and 907 G were first observed by the MIT group using 
a Bose-Einstein condensate of ^^Na atoms [14]. They demonstrated that by 
varying the magnetic field they could change the scattering length by more 
than an order of magnitude. They also observed enhanced inelastic losses near 
the resonances. In Ref. [15], the inelastic losses were studied in greater detail. 

For ^^Rb atoms, the spin-singlet scattering length is large and the spin-triplet 
scattering length is relatively large. They differ from the natural scale £vdw by 
factors of about 17 and —2.4, respectively. Thus most of the hyperfine spin 
states have large scattering lengths at B = 0. However, there are Feshbach 
resonances at which the scattering lengths diverge. For example, the |2, —2) 
hyperfine state has a Feshbach resonance near 155 G. This Feshbach resonance 




(8) 
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was first observed by the Texas group using an ultracold gas of ^^Rb atoms [16]. 
The position and width of this Fcshbach resonance were determined precisely 
by the JILA group [17]. They were used to improve the parameters of the 
atomic potential for Rb and to predict the scattering lengths of the hyperfine 
states of ^^Rb and ^^Rb. They also demonstrated that by varying the magnetic 
field, they could change the collision rate by 4 orders of magnitude and change 
the sign of the scattering length. 



2.4 The resonant and scaling limits 

We have defined a large scattering length to be one that satisfies \a\ S> i, 
where i is the natural low-energy length scale. The corrections to the universal 
behavior are suppressed by powers of i/\a\. There are two obvious limits in 
which the size of these corrections decreases to zero: 

• the resonant or unitary limit a — ±00 with i fixed, 

• the scaling or zero-range limit: £ ^ with a fixed. 

It will sometimes also be useful to consider systems in which the resonant and 
scaling limits are achieved simultaneously: a — ±00 and £ — 0. 

The resonant limit is also sometimes called the unitary limit, because in this 
limit the S-wave contribution to the cross section at low energy saturates 
its unitarity bound cr^^^^^ < Stt/A;^. The resonant hmit can be approached by 
tuning the depth of the interatomic potential to a critical value for which there 
is a 2-body bound state exactly at the 2-body threshold. The resonant limit 
can also be approached by tuning the magnetic field to a Fcshbach resonance. 
Since a — ±00 in the resonant limit, one might expect that in this limit 
the natural low-energy length scale £ is the only important length scale at 
low energies. This is true in the 2-body sector. However, the Efimov effect 
reveals that there is another length scale in the 3-body sector. In the resonant 
limit, there are infinitely many, arbitrarily-shallow 3-body bound states with 
a spectrum of the form 

Et^ — >{e-''^''^'°y'''* h'^Kl/m, asn^+00 with a = ±00, (9) 

where e~'^'^^^° ~ 1/515.03 for identical bosons and be interpreted as the 

approximate binding wave number of the Efimov state labelled by the integer 
n*. If we chose a different integer n*, the value of /t* would change by some 
power of e^/"«. Thus defined by Eq. (9) only up multiplicative factors of 
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The scaling limit is also sometimes called the zero-range limit, because it can 
be reached by simultaneously tuning the range of the potential to zero and its 
depth to oo in such a way that the scattering length is fixed. The terminology 
"scaling hmit" seems to have been first used in Ref . [18] . The scahng limit may 
at first seem a little contrived, but it has proved to be a powerful concept. 
It can be defined by specifying the phase shifts for 2-body scattering. In the 
scaling limit, the S-wave phase shift So{k) has the simple form 



kcotSoik) ^ -1/a, (10) 

and the phase shifts for all higher partial waves vanish. In the scaling 

hmit, the scattering length a sets the scale for most low-energy observables. It 
is the only length scale in the 2-body sector. However, as we shall see, in the 
3-body sector, observables can also have logarithmic dependence on a second 
scale. In the scaling limit, there are infinitely many arbitrarily-deep 3-body 
bound states with a spectrum of the form [19,20] 



4") ^ Ce-2./.or-"* ^ asn^-oc with £^0. (11) 

Thus the spectrum is characterized by a parameter with dimensions of wave 
number. 

The scaling limit may appear to be pathological, because the spectrum of 3- 
body bound states in Eq. (11) is unbounded from below. However, the deep 3- 
body bound states have a negligible effect on the low-energy physics of interest. 
The pathologies of the scaling limit can be avoided simply by keeping in mind 
that the original physical problem before taking the scaling limit had a natural 
low-energy length scale £. Associated with this length scale is a natural energy 
scale jmi^. Any predictions involving energies comparable to or larger than 
the natural low-energy length scale are artifacts of the scaling limit. Thus 
when we use the scaling limit to describe a physical system, any predictions 
involving energies \E\ > fi' lm£? should be ignored. 

In spite of its pathologies, we shall take the scaling limit as a starting point 
for describing atoms with large scattering length. We will treat the devia- 
tions from the scaling limit as perturbations. Our motivation is that when the 
scattering length is large, there are intricate correlations between 3-body ob- 
servables associated with the Efimov effect that can be easily lost by numerical 
approximations. By taking the scaling limit, we can build in these intricate 
correlations exactly at high energy. Although these correlations are unphysical 
at high energy, this does not prevent us from describing low-energy physics 
accurately. It does, however, guarantee that the intricate 3-body correlations 
are recovered automatically in the resonant limit a — > ±oo. 
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2.5 Hyperspherical coordinates 



The universal aspects of the 3-body problem can be understood most easily 
by formulating it in terms of hyperspherical coordinates. A good introduction 
to hyperspherical coordinates and a thorough review of the hyperspherical 
formalism is given in a recent review by Nielsen, Fedorov, Jensen, and Garrido 
[21]. 

In order to define hyperspherical coordinates, we first introduce Jacobi co- 
ordinates. A set of Jacobi coordinates consists of the separation vector rij 
between a pair of atoms and the separation vector Vk^ij of the third atom from 
the center-of-mass of the pair. For atoms of equal mass, the Jacobi coordinates 
are 

rij = T-i - Vj , rk,ij = r-fc - + Tj- ) . (12) 

The hyperradius R is the root-mean-square separation of the three atoms: 

= I (rj, + rl, + rl) = irj + |r^^,^. . (13) 

The hyperradius is small only if all three atoms are close together. It is large 
if any single atom is far from the other two. The Delves hyperangle [22] ak is 
defined by 



ak = arctan ( ^^'^ \ , (14) 

where {i,j, k) is a permutation of (1, 2, 3). The range of the hyperangle ak is 
from to iyr. It is near when atom k is far from atoms i and j, and it is 
near ^tt when atom k is near the center of mass of atoms i and j. 

The Schrodinger equation for the stationary wave function ^(t*i, t"2, ra) of 
three atoms with mass m interacting through a potential V is 



E V- + V{n, r2, ra) j ^ = E^> . (15) 

The wave function ^ in the center-of-mass frame depends on 6 independent 
coordinates. A convenient choice consists of the hyperradius i?, one of the 

hyperangles a^, and the unit vectors fjj and Vksj- We will refer to the 5 
dimensionless variables {ak,'f'ij,rk,ij) as hyperangular variables and denote 
them collectively by Q. When expressed in terms of hyperspherical coordinates, 
the Schrodinger equation reduces to 
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(16) 



where is the hyperradial kinetic energy operator, 



^ 2m 



5 d 



RdR 



2m 
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(17) 



and Th is the kinetic energy operator associated with the hyperangular vari- 
ables. 

The Faddeev equations are a set of three differential equations equivalent to 
the 3-body Schrodinger equation that exploit the simplifications associated 
with configurations consisting of a 2-body cluster that is well-separated from 
the third atom. The solutions to the Faddeev equations are three Faddeev 
wavefunctions whose sum is a solution to the Schrodinger equation: 



*(r-i, T2, r-a) = V^^H^23, ^,23) + V'^^H^si, ^"2,31) + '^^^^^2, ^3,12) ■ 



(18) 



We restrict our attention to states with total angular momentum quantum 
number L = 0. At low energies, we can make an additional simplifying as- 
sumption of neglecting subsystem angular momenta. If the three particles are 
identical bosons, the three Faddeev wave functions can be expressed in terms 
of a single function ip{R, a): 



*(n, r-2, T-a) = TpiR, ai) + ip{R, cva) + tp{R, as) . 



(19) 



The three Faddeev equations can be reduced to a single integro-differential 
equation for the Faddeev wavefunction ■?/'(-R, «)■ (See Ref. [6] for more details.) 

A convenient way to solve the resulting equation is to use a hyperspherical 
expansion. For each value of R, the wave function ip{R,a) is expanded in a 
complete set of functions (f)n{R, a) of the hyperangle a: 



ip{R, a) 



i?V2 sin(2a) 



Y,fn{R)MR:C^)- 



(20) 



The functions (f)n{R,a) are solutions to an eigenvalue equation in the hy- 
perangle a that is parametric in the hyperradius R. The eigenvalues Xn{R) 
determine channel potentials for the hyperradial variable: 



Vn{R) = [Xn{R) - 4] 



2mR^ 



(21) 
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The hypcrradial wavcfunctions fn{R) satisfy an infinite set of coupled par- 
tial differential equations. In the adiabatic hyperspherical approximafion [23], 
the coupling terms are neglected and the equations decouple. They reduce to 
independent hyperradial equations for each of the hyperspherical potentials: 



2m 



15 ^ 



(22) 



In the hyperspherical close- coupling approximation [23], the diagonal coupling 
terms arc also retained. This approximation is more accurate, because it is 
variational in character. 



3 Two Identical Bosons 



In this section, we summarize the universal properties of two identical bosons 
in the scaling hmit in which the large scattering length is the only length scale. 



3.1 Atom-atom scattering 



The cross section for low-energy atom-atom scattering is a universal function 
of the scattering length and the collision energy E. By low energy, we mean 
E — h'^k'^/m much smaller than the natural low-energy scale h'^/mP. The 
partial wave expansion in Eq. (4) expresses the scattering amplitude in terms 
of phase shifts Siik). In the scaling limit, all the phase shifts vanish except for 
So{k) and all the coefficients in the low-energy expansion of kcotSo{k) vanish 
with the exception of the leading term — 1/a. The scattering amplitude in 
Eq. (4) reduces to 



AW = =T7J3«. (23) 



and the differential cross section in Eq. (3) is 



1 + ■ ^ > 

The cross section is obtained by integrating over the solid angle 27r. For wave 
numbers k ^ l/|a|, the differential cross section has the scale-invariant form 
4//^^, which saturates the upper bound from partial-wave unitarity in the 
L = channel. 
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The wave function for atom-atom scattering states at long distances is a uni- 
versal function of the scattering length and the separation r. By long distances, 
we mean r much larger than the natural low-energy length scale £. The sta- 
tionary wave function in the center-of-mass frame for two atoms in an L = 
state with energy E — h'^k'^/m is 



i'AA{f) — - sin [kr — arctan(/ca)] . (25) 
r 

This wave function satisfies the boundary condition 



V'aaW ^C(---) , asr^O, (26) 



r a 



where C is an arbitrary normalization constant. In the scaling limit i ^ 0, the 
effects of short-distances r ~ £ enter only through this boundary condition. 



3.2 The shallow dimer 



The spectrum of shallow 2-body bound states is also universal. By a shallow 
bound state, we mean one whose binding energy Ed is much smaller than the 
natural low-energy scale h'^/mi'^. For a < 0, there are no shallow bound states. 
For a > 0, there is a single shallow bound state, which we will refer to as the 
shallow dimer, or simply as the dimer for brevity. The binding energy Ed of 
the dimer in the scaling limit is 



E 



D 



(27) 



The wave function of the dimer at separations r ^ £ is also universal: 



V'D(r-) = -e-'-ZV (28) 
r 

The size of the dimer is roughly a. A quantitative measure of the size is the 
mean separation of the atoms: (r) = a/2. The wavefunction in Eq. (28) satisfies 
the boundary condition in Eq. (26). 
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3.3 Continuous scaling symmetry 



The only parameter in the universal expressions for the cross-section in Eq. (24) 
and the binding energy in Eq. (27) is the scattering length a. The fact that 
low-energy observables depend only on a single dimensionful parameter a can 
be expressed formally in terms of a continuous scaling symmetry that consists 
of rescaUng a, the coordinate r, and the time t by appropriate powers of a 
positive number A: 

a — >\a, r — > \r , t — > \H . (29) 

The scaling of the time by the square of the scaling factor for lengths is 
natural in a nonrelativistic system. Under this symmetry, observables, such as 
the dimer binding energy or the atom-atom cross section gaai scale with 
the powers of A suggested by dimensional analysis. 

The scaling symmetry strongly constrains the dependence of the observables 
on the scattering length and on kinematic variables. As a simple example, 
the dimer binding energy scales as Ed — > \~'^Ed- The scaling symmetry con- 
strains its dependence on the scattering length to be proportional to 
in agreement with the explicit formula in Eq. (27). As another example, the 
atom-atom cross section scales as aAA ^cfaa- The scaling symmetry con- 
strains its dependence on the energy and the scattering length: 

aAA{X-^E; Xa) = X^aAA^E- a) . (30) 

The explicit expression for the differential cross section in Fig. (24) is consis- 
tent with this constraint. 

The scattering length a changes discontinuously between -|-oo and — oo as the 
system is tuned through its critical point. Since 1/a changes smoothly, this is 
a more convenient interaction variable. To exhibit the scaling symmetry most 
clearly, it is convenient to use an energy variable with the same dimensions as 
the interaction variable. A convenient choice is the wave number variable 

K = sign{E){m\E\/n^f''^. (31) 

The set of all possible low-energy 2-body states in the scaling limit can be 
represented as points {ar^ , K) on the plane whose horizontal axis is 1/a and 
whose vertical axis is K. It is convenient to also introduce polar coordinates 
consisting of a radial variable H and an angular variable ^ defined by 

l/a^Hcosi, K^Hsm^. (32) 
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In terms of these polar coordinates, the scahng symmetry given by Eqs. (29) 
is simply a rescaling of the radial variable: H X~^H. 



K 



n 



AA 



AA 




^^^^^^^ 1/a 



Fig. 1. The a -K plane for the 2-body problem. The allowed region for atom-atom 
scattering states are the two quadrants labelled AA. The heavy line labelled D is 
the shallow dimer. The cross-hatching indicates the 2-atom threshold. 

The a~^-K plane for the 2-body system in the scaling limit is shown in Fig. 1. 
The possible states are atom-atom scattering states {AA) and the shallow 
dimer [D). The threshold for atom-atom scattering states is indicated by the 
hatched area. The shallow dimer is represented by the heavy line given by the 
ray ^ = — ^tt. A given physical system has a specific value of the scattering 
length, and so is represented by a vertical line, such as the dashed line in 
Fig. 1. Changing a corresponds to sweeping the line horizontally across the 
page. The resonant limit corresponds to tuning the vertical line to the K axis. 

3.4. Scaling violations 

The continuous scaling symmetry is a trivial consequence of the fact that a is 
the only length scale that remains nonzero in the scaling limit. For real atoms, 
the scaling limit can only be an approximation. There arc scaling violations 
that give corrections that are suppressed by powers oi £/\a\. In the 2-body 
sector, the most important scaling violations come from the S-wave effective 
range defined by the effective-range expansion in Eq. (5). 

The differential cross section for atom-atom scattering can be expanded in 
powers oiija with ok fixed: 
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da Ao? 



(33) 



dn 1 + a?k'^ 



1 + - 



a 1 + a?k'^ 



+ ... 



The leading term is tiie universal expression in Eq. (24). The ncxt-to-leading 
term is determined by Ts. 

If a > 0, we can also consider the scaling violations to the binding energy of 
the shallow dimer. The binding energy can be expanded in powers of £/a: 



The leading term is the universal expression in Eq. (27). The first two correc- 
tion terms are determined by r^. 



4 Three Identical Bosons 

In this section, we summarize the universal properties of three identical bosons 
in the scaling limit in which the only scales are those provided by the large 
scattering length a and the Efimov parameter 



4.1 Boundary condition at short distances 

When the scattering length is large compared to the natural low-energy length 
scale ^, the range of the hyperradius it! includes four important regions. It is 
useful to give names to each of these regions: 

• the short- distance region i? < |£|, 

• the scale-invariant region £ <C -R <^ |a|, 

• the long-distance region ~ |a|, 

• the asymptotic region \a\. 

In the scaling limit £ — > 0, the short-distance region shrinks to 0. Its effects 
can however be taken into account through a boundary condition on the hy- 
perradial wavefunction in the scale-invariant region. 

The hyperspherical expansion of the Faddeev wave function is given in Eq. (20). 
In the scaling limit, the channel eigenvalues A„(i?) determined by the eigen- 
value equation for the hyperangular functions (pniR, ct) satisfy [2] 




(34) 
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COS (A^/^f ) - sin [X^^) = sin [x'/^) ^ . (35) 

The lowest channel potentials Vn{R) defined by Eq. (21) are shown in Fig. 2. 
For R <^ \a\, the channel eigenvalues approach constants A„(0) and the channel 




Fig. 2. The three lowest hyperspherical potentials Vn{R) scaled by h^/ma'^ for a > 
(solid lines) and for a < (dashed lines). 

potentials in Eq. (21) have the scale-invariant behavior 1/R^. Whether a is 
positive or negative, the lowest eigenvalue Ao(-R) is negative at i? = 0. It 
can be expressed as Ao(0) = — Sq, where sq ~ 1.00624 is the solution to the 
transcendental equation 



sn cosh — ^ = sinh . (36) 
2 73 6 ^ ' 

Thus the potential Vo(-R) is an attractive 1/i?^ potential for R <^ \a\: 



Vo{R)^-{4 + sl)^^^, £«i?«|a|. (37) 

All the channel eigenvalues for n > 1 are positive at i? = 0: A„(0) > 19.94. 
Thus all the channel potentials Vn{R) for n > 1 are repulsive l/R"^ potentials in 
the region R <^ \a\. The hyperradial wave functions fn{R) for ^ > 1 therefore 
decrease exponentially for i? <C |a|. 

In the scaling limit £ — > 0, the lowest adiabatic hypersherical potential in 
Eq. (37) behaves like l/R"^ all the way down to i? = 0. Such a potential 
is too singular to have well-behaved solutions. If the exact solution f{R) 
for the hyperradial wave function at short distances was known, it could be 
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matched onto the solution in the scaling Umit by choosing a hyperradius Rq 
in the scale-invariant region and demanding that the logarithmic derivatives 
Rof'{Ro)/f{R(yj match at that point. If we also choose Rq -C {m\E\/h'^)~^/'^, 
the energy eigenvalue E in Eq. (22) can be neglected relative to the channel 
potential. The most general solution for the hyperradial wave function in the 
scale-invariant region is 



where A and B are arbitrary coefficients. To make the arguments of the loga- 
rithms dimensionless, we have inserted factors of H, the wavenumber variable 
defined in Eq. (32). The terms with the coefficients A and B represent an 
outgoing hyperradial wave and an incoming hyperradial wave, respectively. If 
\A\ < \B\, there is a net flow of probability into the short-distance region. As 
will be discussed in detail in Section 5, such a flow of probability is possible if 
there are deep diatomic molecules. In this section, we assume that there are no 
deep 2-body bound states. The probability in the incoming hyperradial wave 
must therefore be totally reflected at short distances. This requires \A\ — \B\, 
which implies that A and B differ only by a phase: 



A = -e^i^*B (39) 

for some angle 9^. This angle can be expressed as = — Sq ln(if/Ao), where 
Ao is the product of l/i?o and a comphcated function of Rof'{Ro)/ f{Ro). The 
effects of the short- distance region on f{R) enter only through the wavenum- 
ber variable Aq. By solving the hyperradial equation for the scale-invariant 
potential in Eq. (37), we flnd that Aq differs from the 3-body parameter 
defined by the Efimov spectrum in the resonant limit only by a multiplicative 
numerical constant. Thus the angle 9^ in Eq. (39) can be expressed as 



9^^ -Soln{cH/K,) , (40) 
where H is defined in Eq. (32) and c is a constant. 



4-2 Discrete scaling symmetry 



The boundary condition in Eq. (39) gives logarithmic scaling violations that 
give corrections to the scaling limit that are functions of ln(if/«;*). Logarith- 
mic scaling violations do not become less important as one approaches the 
scaling limit, and therefore cannot be treated as perturbations. Because the 
boundary condition enters through the phase in Eq. (39), the logarithmic scal- 
ing violations must be log-periodic functions of H/k^ with period tt/sq- 
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The 3-body sector in the scahng hmit has a trivial continuous scahng symme- 
try defined by Eqs. (29) together with A~^k^,. However, because of the 
log-periodic form of the logarithmic scaling violations, it also has a nontriv- 
ial discrete scaling symmetry. There is a discrete subgroup of the continuous 
scaling symmetry that remains an exact symmetry in the scaling limit: 

— a — ^AqO, r — > XqV , t — > Ag^t , (41) 

where n is an integer, Aq = e'^/*°, and Sq ~ 1.00624 is the solution to the tran- 
scendental equation in Eq. (36). The numerical value of the discrete scaling 
factor is e'^/^° ~ 22.7. Under the discrete scaling symmetry, 3-body observ- 
ables, such as binding energies and cross sections, scale with the integer powers 
of Aq suggested by dimensional analysis. By combining the trivial continuous 
scaling symmetry with the discrete scaling symmetry given by Eqs. (41), we 
can see that k^, is only defined modulo multiplicative factors of Aq. 

The discrete scaling symmetry strongly constrains the dependence of the ob- 
servables on the parameters a and and on kinematic variables. For example, 
the scaling of the atom-dimer cross section is aAD — ^ (A^)^cr^_D- The discrete 
scaling symmetry constrains its dependence on a, n^,, and the energy E: 

aAoiXf'^E- A™a, k,) = X^^aAD^E- a, k,) , (42) 

for all integers m. hi E = 0, the cross section is simply aAD = 47r|aADp, 
where qad is the atom-dimer scattering length. The constraint in Eq. (42) 
implies that the atom-dimer scattering length aAD is proportional to a with a 
coefficient that is a log-periodic function of an^ with period tt/sq. The explicit 
expression for the atom-dimer scattering length is given in Eq. (61) and it is 
indeed consistent with this constraint. 

To illustrate the discrete scaling symmetry, it is convenient to use the inter- 
action variable 1/a and the energy variable K defined in Eq. (31). For a given 
value of K^,, the set of all possible low-energy 3-body states in the scahng limit 
can be represented as points (a~^, K) on the plane whose horizontal axis is 1/a 
and whose vertical axis is K. The discrete scaling transformation in Eqs. (41) 
is simply a rescaling of the radial variable H defined in Eq. (32) with and 
e fixed: H Xq'^H. 

The a~^-K plane for three identical bosons in the scaling limit is shown in 
Fig. 3. The possible states are 3-atom scattering states (AAA), atom-dimer 
scattering states {AD), and Efimov trimers (T). The threshold for scattering 
states is indicated by the hatched area. The Efimov trimers are represented by 
the heavy lines below the threshold!^ There are infinitely many branches of 

^ The curves for the trimer binding energies in Fig. 3 actually correspond to plotting 
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Fig. 3. The a -K plane for the 3-body problem. The allowed regions for 3-atom 
scattering states and for atom-dimer scattering states are labelled AAA and AD, 
respectively. The heavy lines labeled T are three of the infinitely many branches of 
Efimov states. The cross-hatching indicates the threshold for scattering states. 

Efimov trimers, but only a few arc shown. They intercept the vertical axis at 
the points K — — (e'~'^/*°)"~"*K*. A given physical system has a specific value 
of the scattering length, and so is represented by a vertical fine. The resonant 
limit corresponds to tuning the vertical line to the K axis. 

Changing 1/a continuously from a large positive value to a large negative 
value corresponds to sweeping the vertical dashed line in Fig. 3 from right 
to left across the page. The Efimov trimers appear one by one at the atom- 
dimer threshold at positive critical values of 1/a that differ by powers of 
g7r/so ~ 22.7 until there are infinitely many at 1/a = 0. As 1/a continues to 
decrease through negative values, the Efimov trimers disappear one by one 
through the 3-atom threshold at negative critical values of 1/a that differ 
by powers of e"^/*". We will focus on the specific branch of Efimov trimers 
labelled by the integer n = n*, which is illustrated in Fig. 4. At some positive 
critical value a = a*, this branch of Efimov trimers appears at the atom-dimer 
threshold: E^^*^ = Ed- As a increases, its binding energy Ej^*^ — E^) relative 
to the atom-dimer threshold increases but its binding energy E^*^ relative to 
the 3-atom threshold decreases monotonically. As a — > oo, the binding energy 
approaches a nontrivial limit: e'^*'' — > h^K,l/m. For a < 0, as \a\ decreases. 



H^/^sinS^ versus i^^/'^cos^. This effectively reduces the discrete scaling factor 22.7 
down to 22.7^/^ = 2.2, allowing a greater range of and K to be shown in the 
Figure. 



21 





'(3 J); 




< 


► 




1/a 



Fig. 4. The energy variable K for the branch of Efimov trimers labelled by n = 
as a function of 1/a. In the resonant limit, the binding wave number is The 
branch disappears through the atom-dimer threshold at a = a* and through the 

3- atom threshold at a = a'^. 

e!P*^ continues to decrease monotonically. Finally, at some negative critical 
value a — a'^, it disappears through the 3-atom threshold. 

4- 3 Efimov trimers 

The binding energies of the Efimov states are functions of a and k^,. Efimov 
showed that the calculation of the binding energies for all the Efimov states 
could be reduced to the calculation of a single universal function of ^. The 
Efimov states can be interpreted as bound states in the lowest adiabatic hy- 
perspherical potential. This potential has a scale-invariant region where the 
general solution is the sum of an outgoing hyperradial wave and an incom- 
ing hyperradial wave as in Eq. (38). Bound states occur at energies for which 
the wave reflected from the long-distance region i? ~ |a| come into resonance 
with the wave reflected from the short-distance region i? ~ The resonance 
condition can be expressed as 

26^ + A(^) = mod 27r , (43) 

where A(,^)/2 is the phase shift of a hyperradial wave that is reflected from the 
long-distance region. Using the expression for in Eq. (40) and the deflnitions 
for H and ^ in Eq. (32), we obtain Efimov's equation for the binding energies: 
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+ ^ = (e-^-/-)"-'^' exp [A(0/.o] ^ , (44) 
where the angle ^ is defined by 

tan^ = -{mET/n'^y^^ a. (45) 

We have absorbed the constant c in Eq. (40) into the function A(^) so that it 
satisfies A(— ^tt) — 0. 

Once the universal function A(^) has been calculated, the binding energies for 
all the Efimov states for any values of a and be obtained by solving 

Eq. (44). The equation is the same for different Efimov states except for the 
factor of (e~27r/sojn right side. 

Efimov's universal function A(^) was calculated in Ref. [24] with a few digits 
of precision over the entire range — tt < ^ < It has been calculated in 

Ref. [25] with a precision of about 12 digits for a > 0, which corresponds to 
the range — ^tt < ^ < — |vr. It decreases from 6.0273 at ,^ = — ^tt to at 
^ = — ^TT and then to about —0.89 at ^ = — tt. A parameterization of A(^) 
is given in Ref. [6] that has errors less than about 0.01, at least in the range 
-|7r<e< -|7r. 

In the resonant limit a — ±oo, the spectrum of the Efimov states is 
particularly simple. In this limit, ^ — > —^n and A(^) — > 0, so the solutions to 
Eq. (44) reduce to 



4n) ^ / -2./.or- ^ a = ±oo, (46) 

where is the binding wave number for the Efimov state labeled by n = n*. 
The spectrum in Eq. (46) is geometric, with the binding energies of successive 
Efimov states having the ratio e~'^'^l^° 1/515.03. The Schrodinger wave 
function in the center-of-mass frame for an Efimov state with binding energy 
Et = h^K^/m is 



3 sinh snff — ttj 

*(n, r„ ra) = R-"'h{R) E 

7^1 sm(2ai) 



(47) 



The hyperradial wavefunction is 



fo{R) = R'/^Ki,,{V2KR), 



(48) 
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where KisQ {z) is a Bessel function with imaginary index. A quantitative mea- 
sure of the size of a 3-body bound state is the mean-square hyperradius: 

= (e^-/-)"-"* . (49) 

Thus the root-mean-square hyperradius for each successively shallower Efimov 
state is larger than the previous one by e^^^° a; 22.7. 

The Efimov trimers disappear through the atom-dimer threshold at positive 
critical values of a, as illustrated in Fig. 4. For the branch of Efimov trimers 
labelled by n = n*, the critical value is 



a* = 0.0707645086901 . (50) 

The other critical values are (e'^/*°)"a*, where n is an integer. The binding 
energy Ej^^ for the Efimov state just below the atom-dimer threshold is 



E^T^ ^ED\i + 0.164 W{a/a,)\ . (51) 
The errors in this approximation scale as (a — a*)^. 

Efimov states near the atom-dimer threshold can be understood intuitively as 

2- body systems composed of an atom of mass m and a dimer of mass 2m. We 
can exploit the universality of the 2-body systems with large scattering lengths 
to deduce some properties of the shallowest Efimov state when it is close to 
the atom-dimer threshold. The analog of the universal formula in Eq. (27) is 
obtained by replacing the reduced mass m/2 of the atoms by the reduced mass 
2m/3 of the atom and dimer and by repacing a by the atom-dimer scattering 
length a AD which diverges at a = a*. Thus the binding energy relative to the 

3- atom threshold can be approximated by 



E^^^ ^Ed + . (52) 

4ma, 



An explicit expression for the atom-dimer scattering length qad is given in 
Eq. (61). The binding energy in Eq. (52) agrees with the result in Eq. (51) up to 
errors of order a/a\j^, which scales as (a — a*)^. We can also use universality to 
deduce the wave function of the Efimov trimer. The Schrodinger wave function 
can be expressed as the sum of three Faddeev wave functions as in Eq. (19). 
In the limit ri 23 ^ r23, the first Faddeev wave function should have the form 



'0^^n^23, ^"1,23) ~ ^A£)(ri,23) i>D{r23) , (53) 
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where ipDif) is the dimer wave function given in Eq. (28) and iljAoi'T') is the 
analogous universal wave function for a shallow bound state consisting of two 
particles with large positive scattering length a ad'- 



ipAoir) = -e-'-/'^^- . (54) 
r 

This Faddeev wave function can be expressed in terms of hyperspherical co- 
ordinates using Eqs. (13) and (14). Most of the support of the probability 
density j^'p is concentrated in the region in which the hyperradius is very 
large, R ~ a ad, and one of the three hyperangles is very small, CKj <^ 1. The 
mean-square hyperradius can be calculated easily when qad ^ o,: 



This result can be obtained more easily simply by using the universal atom- 
dimer wave function in Eq. (54) and the approximate expression R^ for 
the hyperradius. 

The Efimov trimers disappear through the 3-atom threshold at negative crit- 
ical values of a, as illustrated in Fig. 4. For the branch of Efimov trimers 
labelled hy n — n^, the critical value is 



< = -1.56(5) (56) 

The other critical values are (e'^/*")"a'^, where n is an integer. In contrast to 
a* in Eq. (50), only a few of digits of precision are currently available for a'^. 
Comparing Eqs. (50) and (56), we see that a'^ ^ —22.0 a*. 

There is an Efimov trimer at the 3-atom threshold when a has the negative 
critical value a'^. As can be seen in Fig. 4, the binding energy e!^^ of the Efimov 
trimer increases rapidly as a function of a when \a\ exceeds \a'^\. The increase 
is so rapid that in Ref. [24], a parameterization of the universal function A(^) 
in Eq. (44) with an essential singularity at a = a'^ was used to get a good fit 
to the binding energy. The result of Refs. [26] and [27] seem to indicate that 
E^^ is actually linear in 1/a near a = a'^ with a large slope. The results of a 
calculation in Ref. [27] can be used to obtain the approximation 



E, 



(n) 



1.1 In 



(57) 



When \a\ decreases below \a'^\, the Efimov trimer does not immediately dis- 
appear, but instead becomes a resonance that decays into three atoms. This 
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resonance is associated with a pole at a complex energy Ej-es — i^ies/^ with a 
positive real part and a negative imaginary part. This complex energy is the 
analytic continuation of the energy —Ejf^ of the shallowest Efimov trimer for 
\a\ < \a'J\. For a close to a'^, the real part of the resonance energy is given 
simply by the negative of the expression in Eq. (57). The authors of Ref. [27] 
gave a parameterization of the imaginary part of the resonance energy that 
scales as (a — a'^Y near the threshold. However their numerical results seem 
to suggest that Fres scales as a higher power of a — a'^. 

There is a common misconception in the literature that Efimov states must 
have binding energies that differ by multiplicative factors of 515.03. However, 
this ratio applies only in the resonant limit a ±oo. The ratio 4" ^V-^t"'' 
of the binding energies of adjacent Efimov trimers can be much smaller than 
515 if a > and much larger than 515 if a < 0. The smallest ratios occur 
at the critical values (e'^/*°)"a*, where a* is given in Eq. (50). The accurate 
results of Ref. [25] for the binding energies Et of the first few Efimov states 
in units of E^ — h^/ma^ are 



E^'> ^Ed, (58a) 
4^-1) = 6.75029015026 Ed , (58b) 
4^-2) = 1406. 13039320 Ed- (58c) 

Thus, if a > 0, the ratio E^~^^ / E^'' of the binding energies for the two 

shallowest Efimov trimers can range from about 6.75 to about 208. The largest 
ratios occur at the critical values (e'^/*'")"a'^, where a'^ is given in Eq. (56). The 
binding energies E^ of the first few Efimov states are [24] 



E^T> =0, (59a) 
4^-1) = 1.09 X 10=^ h'^/ma'^ , (59b) 
= 5.97 X 10^ h^/ma'^ . (59c) 

Thus, if a < 0, the ratio Ej^~^'' /Ej^^ of the binding energies for the two 
shallowest Efimov states can range from about 550 to oo. 



4-4 Atom-dimer elastic scattering 



The differential cross section for elastic atom-dimer scattering near the atom- 
dimer threshold can be expressed in terms of the atom-dimer scattering length 
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da AD 



as E ^ —Ed ■ 



(60) 



The discrete scaling symmetry implies that a ad jo. must be a log-periodic 
function of a/t* with period tt/sq- Its functional form was deduced by Efimov 
up to a few numerical constants. The numerical constants were first calculated 
in Ref. [28]. The atom-dimer scattering length is 



where a* ~ 0.071 k"^ is given to high precision in Eq. (50). The atom-dimer 
scattering length diverges if a has one of the values (e'^/*°)"'a* for which there 
is an Efimov state at the atom-dimer threshold. It vanishes if a has one of the 
values (e^/"o)"0.38a*. 

^.5 Three-body recombination 

Three-body recombination is a process in which three atoms coUide to form a 
diatomic molecule and an atom. The energy released by the binding energy of 
the molecule goes into the kinetic energies of the molecule and the recoiling 
atom. The 3-body recombination rate depends on the momenta of the three 
incoming atoms. If their momenta are sufficiently small compared to l/|a|, 
the dependence on the momenta can be neglected, and the recombination 
rate reduces to a constant. The recombination event rate constant a is defined 
such that the number of recombination events per time and per volume in a 
gas of cold atoms with number density ua is an\. If the atom and the dimer 
produced by the recombination process have large enough kinetic energies to 
escape from the system, the rate of decrease in the number density of atoms 
is 

= -San^ , (62) 

In a Bose-Einstein condensate, the three atoms must all be in the same quan- 
tum state, so the coefficient of n\ in Eq. (62) must be multiplied by 1/3! 
to account for the symmetrization of the wave functions of the three iden- 
tical particles [29]. This prediction was first tested by the JILA group [30]. 
They measured the 3-body loss rates in an ultracold gas of ^^Rb atoms in the 
If.rrif) = |1, —1) hyperfine state, both above and below the critical tempera- 
ture for Bose-Einstein condensation. They found that the loss rate was smaller 
in the Bose-Einstein condensate by a factor of 7.4(2.6), in agreement with the 
predicted value of 6 [29] . 




(61) 



27 



If the scattering length a is negative, the molecule can only be a deep di- 
atomic molecule with binding energy of order h^/mP or larger. However, if 
a is positive and unnaturally large (a ^ i), the molecule can also be the 
shallow dimer with binding energy Ejj = /mo?. Three-body recombination 
into deep dimers will be discussed in Section 5. In this section, we assume 
there arc no deep dimers. We therefore assume a > and focTis on 3-body 
recombination into the shallow dimer. We denote the contribution to the rate 
constant a from 3-body recombination into the shallow dimer by ttshaiiow 

Dimensional analysis together with the discrete scaling symmetry implies that 
^shallow is proportional to ha'^/m with a coefficient that is a log-periodic func- 
tion of a/t* with period tt/sq. An analytic expression for ashaiiow has recently 
been derived [31,32]: 



1287r2(47r-3V3) sin^[so ln(a/a^,o)] ha"^ , , 

•^shallow — . ,2/ \~ 9? 1 7 '/ vi ' V""^J 

smh (ttsoJ -l- cos^[so ln(a/a*ojJ ^ 

where a*o differs from by a multiplicative constant that is known only to 
a couple of digits of accuracy: 



a*o ~ 4.5 a* 0.32 (64) 
The maximum value of the coefficient of ha^ /m in Eq. (63) is 



_ 1287r^(47r - 3^3) 

smh (TTSq) 

Its numerical value is Cmax = 67.1177. We can exploit the fact that sinh^(7rso) ~ 
139 is large to simplify the expression in Eq. (63). The rate constant can be 
approximated with an error of less than 1% of Cmax.^a'^ / m by 



dishallow ~ 67.12 sin^[so ln(a/a*o)] fia^/m . (66) 

This approximate functional form of the rate constant was first deduced in 
Refs. [33,34]. The coefficient Cmax and the relation between a*o and was 
calculated accurately in Refs. [35,36]. 

The most remarkable feature of the analytic expression in Eq. (63) and the 
approximate expression in Eq. (66) is that the coefficient of ha'^/m oscillates 
between zero and 67.12 as a function of ln(a). In particular, ^shallow has zeroes 
at values of a given by (e'^/*")" a^oi where a^,o ~ 4.5 a*. The maxima of ttshaiiow 
in Eq. (66) occur at the values of a for which tan[so ln(a/a*o)] = —Sq/2, which 
are a ^ (e'^/'^o)" 64.3 a* ^ (e''/*o)M.6 fi;-^ 
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Since the zeroes in ^shallow are so remarkable, it is worth enumerating the 
effects that will tend to fill in the zeroes, turning them into local minima 
of ashaiiow The zeroes arise from the interference between two pathways as- 
sociated with the lowest two adiabatic hyperspherical channels. There may 
be additional contributions from coupling to higher hyperspherical channels, 
but in the scaling limit they are strongly suppressed numerically. The zero in 
^shallow for a = 0*0 is exact only at threshold. If the recombining atoms have 
collision energy E, ctshaiiow goes to zero as as E [37]. Thus thermal 
effects that give nonzero collision energy to the atoms will tend to fill in the 
zeroes. Finally, if the 2-body potential supports deep diatomic molecules, their 
effects will tend to fill in the zeros of ashaiiow as described in Section 5.4. Fur- 
thermore, as described in Section 5.4, 3-body recombination into those deep 
dimers gives an additional nonzero contribution a^eep to the rate constant a. 



5 Effects of Deep Diatomic Molecules 

In this Section, we describe the effects of deep diatomic molecules on the 
universal properties for three identical bosons in the scaling limit. 

5.1 Boundary condition at short distances 

The existence of deep dimers requires a modification of the boundary condi- 
tion on the hyperradial wave function at short distances. The general solution 
for the hyperradial wave function in the lowest hyperspherical potential in the 
scale- invariant region £ -C -C \a\ is given in Eq. (38). The boundary condi- 
tion in Eq. (39), which takes into account the effects of short distances -R < 
follows from the assumption that all the probability in an incoming hyper- 
radial wave is refiected back from the short-distance region into an outgoing 
hyperradial wave. If there are deep dimers, some of the probability in the in- 
coming hyperradial wave that fiows into the short-distance region emerges in 
the form of scattering states that consist of an atom and a deep dimer with 
large kinetic energy but small total energy. We will refer to these states as 
high-energy atom- dimer scattering states. 

The 2-body potentials for the alkali atoms other than hydrogen support di- 
atomic molecules with many vibrational energy levels. If it was necessary to 
take into account each of these deep dimers explicitly, the problem would be 
extremely difficult. Fortunately the cumulative effect of all the deep dimers 
on low-energy 3-atom observables can be taken into account through one ad- 
ditional parameter: an inelasticity parameter 7]^ that determines the widths 
of the Efimov trimers. In the scaling limit, the low-energy 3-body observables 
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are completely determined by a, k*, and 77*. 

The reason the cumulative effects of the deep dimers can be described by a 
single number 77^, is that all pathways from a low-energy 3-atom state to a 
high-energy atom-dimer scattering state must flow through the lowest hyper- 
spherical potential, which in the scale-invariant region has the form given in 
Eq. (37). The reason for this is that in order to reach a high-energy atom-dimer 
scattering state, the system must pass through an intermediate configuration 
in which all three atoms are simultaneously close together with a hyperradius 
R of order i or smaller. Such small values of R are accessible to a low-energy 
3-atom state only through the lowest hyperspherical potential. 

If there are deep dimers, some of the probability in a hyperradial wave that 
flows to short distances emerges in the form of high-energy atom-dimer scat- 
tering states. We denote the fraction of the probability that is reflected back 
to long distances by e~^^* . We refer to t/^, as the Efimov width parameter. The 
amplitude A of the outgoing wave in Eq. (38) then differs from the ampli- 
tude B of the incoming wave not only by a phase as in Eq. (39) but also by 
a suppression factor e"^''*. Thus if there are deep 2-body bound states, the 
boundary condition in Eq. (39) must be replaced by 



A = _e-2'?*+2^e*5. (67) 



The angle is related to the 3-body parameter by Eq. (40). The parameters 
and ?7* appear in the boundary condition in Eq. (67) in the combination 
9^ + ir]^. If the universal results for the case 77* = can be expressed as an 
analytic function of ln(K*), then the corresponding universal results for the 
case in which there are deep 2-body bound states can be obtained simply by 
the substitution ln(«;*) ln(«;*) +ir]^/so- 

One can also obtain universal results for the total probability of transitions 
from low-energy 3-atom or atom-dimer scattering states to high-energy atom- 
dimer scattering states. The transition rate for any individual deep dimer is 
sensitive to the details of wave functions in the short-distance region i? ~ ^. 
However the inclusive rate summed over all deep dimers is much less sensitive 
to short distances, because it is constrained by probability conservation. The 
inclusive rates are given by universal expressions that include a factor of 1 — 
g-4r;, j^ggj^ wMch is the probability that a hyperradial wave is not reflected 
from the short-distance region. 
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5.2 Widths of Efimov states 



One obvious consequence of the existence of deep dimers is that the Efimov 
trimers are no longer sharp states. They are resonances with nonzero widths, 

because they can decay into an atom and a deep dimer. The binding energy Ej' 
and width of an Efimov resonance can be obtained as a complex eigenvalue 
E = —{Et + of the 3-body Schrodinger equation. 

In the absence of deep dimers, the binding energies of Efimov trimers satisfy 
Eq. (44), where A(,^)/2 is the phase shift of a hyperradial wave that is refiected 
from the long-distance region ~ \a\. To obtain the corresponding equation 
in the case of deep dimers, we need only make the substitution ^* — > ^* + irj^ 
in Eq. (43): 



2{0^ + irj^) + A{C) = mod 27r. 



(68) 



This can be satisfied only if we aUow complex values of ^ in the argument of 
A. Using the expression for 6^, in Eq. (40) and inserting the definition of H in 
Eq. (32), we obtain the equation 



2 ma^ ^ ^ 



exp 



A(0 + 2^r7, 
So 



2 ,.2 



m 



(69) 



where the complex- valued angle ^ is defined by 



1 /9 

tan^ = -(m{ET + irT/2)/n^^ a 



(70) 



To solve this equation for Ex and Ft, we need the analytic continuation of 
A(^) to complex values of ^. The parametrization for A(^) in Ref. [6] should 
be accurate for complex values of ^ with sufficiently small imaginary parts, 
except near ^ = — vr where the parametrization used an expansion parameter 
that was not based on any analytic understanding of the behavior near the 
endpoint. If the analytic continuation of A(^) were known, the binding energy 
and width of one Efimov state could be used to determine and 77*. The 
remaining Efimov states and their widths could then be calculated by solving 
Eq. (69). 

If the Efimov width parameter r]^, is extremely small, the right side of Eq. (69) 
can be expanded to first order in 77*. The resulting expression for the width is 




Et + 



ma'' 



(71) 
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For the shallowest Efimov states, the order of magnitude of the width is simply 
rj^.h'^ / ma^ . The widths of the deeper Efimov states are proportional to their 
binding energies, which behave asymptotically like Eq. (9). This geometric 
increase in the widths of deeper Efimov states has been observed in calculations 
of the elastic scattering of atoms with deep dimers [39] . 



5.3 Atom-dimer scattering 



The effects of deep dimers modify the universal expressions for low-energy 
3-body scattering observables derived in Section 4. If the universal expression 
for a scattering amplitude for the case of no deep dimers is expressed as 
an analytic function of ln(«;*), the corresponding universal expression for the 
case in which there are deep dimers can be obtained simply by substituting 
ln(K*) ln{K^) + irj^/sQ. 

The differential cross section for elastic atom-dimer scattering near the atom- 
dimer threshold is still given by Eq. (60), except that the atom-dimer scatter- 
ing length is now complex valued: 

ttAD = (l-46 + 2.15 cot[so ln(a/a*) + ir]^fj a , (72) 
where a* ~ 0.071«;~^ is given in Eq. (50). The elastic cross section reduces to 



(elastic)/ X c . n sin2[so ln(a/a,) + 0.97] + sinh^ rj^ 

(^AD = -^d) = 84.9 M , • V.2 « 

sm [soln(a/a*jJ smh 77* 



(73) 



The coefficient of is shown in Fig. 5 as a function of a for several values of 
?7*. In the limit r]^ ^ 00, the log-periodic dependence on ok* disappears and 
the cross section reduces to 84.9 a^. 

If there are no deep dimers, atom-dimer scattering is completely elastic below 
the dimer-breakup threshold ka = 2/\/3. The existence of deep dimers opens 
up inelastic channels in which an atom and a shallow dimer with low energy 
collide to form an atom and a deep dimer. The large binding energy of the deep 
dimer is released through the large kinetic energies of the recoihng atom and 
dimer. This process is called dimer relaxation. The optical theorem implies 
that the total cross section for inelastic atom-dimer scattering near the atom- 
dimer threshold E = —Ed is proportional to the imaginary part of the atom- 
dimer scattering length in Eq. (72): 



^(inelastic) ^ ^ T ^'^"^ "^^^ " ^^^^ 
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Fig. 5. The elastic cross section for atom-dimer scattering at threshold in units of 
as a function of a/a* for several values of r/*. 

Thus the cross section for dimer relaxation diverges like 1/k as E approaches 
the atom-dimer threshold. 

The event rate (3 for dimer relaxation is defined so that the number of dimer 
relaxation events per time and per volume in a gas of atoms with number 
density ua and dimers with number density is f^nAUo- If the atom and 
the deep dimer produced by the relaxation process have large enough kinetic 
energies to escape from the system, the rate of decrease in the number densities 
is 



—riA = —riD = -PuaUd ■ (75) 
at at 

The event rate f3 can be expressed in terms of a statistical average of the 
inelastic atom-dimer cross section: 



In the low-temperature limit, /3 reduces to 67r^(—Im aAZ))/m, which can be 
written as [38] 



If ?7* is small, the maximum value of (3 occurs when a ~ [1 + sinh^ 77*/ (2sg)]a*. 

The coefficient of ha/m is shown in Fig. 6 as a function of a/a^, for several 
values of rj^. It displays resonant behavior with maxima when the scattering 




(76) 



20.3sinh(2r7,) ha 
sin^[so ln(a/a*)] + sinh i]^ m 



(77) 
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Fig. 6. The dimer relaxation rate constant in units oifiajm as a function of aja^^ 
for several values of r?* . 

length has one of the values (e'^/*°)"a* for which the peak of an Efimov res- 
onance is at the atom-dimer threshold. In the limit ry* — >■ 0, the maximum 
value 40.6 cothr^,, of the coefficient of ha/m diverges. In the limit r^^, — >■ oo, the 
log-periodic dependence of the coefficient on ok* disappears and it reduces to 
the constant 40.6. 



5.4 Three-body recombination 
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Fig. 7. The 3-body recombination rate constants ctshaiiow (with large-amplitude oscil- 
lations) and adeep (with small-amplitude oscillations) in units of tia^/m as functions 
of a/a* for o > and several values of rj^, . 

If there are no deep dimers, the rate constant ctshaiiow for 3-body recombination 
into the shallow dimer has the remarkable form given in Eq. (63), which has 
zeroes at values a that differ by multiples of e^/^° ~ 22.7. If there are deep 
dimers, the analytic expression for the rate constant ashaiiow is 
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_ 1287r^(47r - 3v^)(sm^[go ln(a/a,o)] + sinh^ 77,) ha^ 

Clshallow — -12/ , \ , 9r 1 / / \1 ' \'°) 

smh (ttso + V*) + cos^[so in(a/a*oJJ "t. 

where a*o ~ 4.5 a* is given in Eq. (64). We can use the approximation sinh(7rso+ 
77*) ~ e'^* sinh(7rso) to simphfy the expression in Eq. (78): 



^shallow ~ 67.12 e ^''* (^sin^[soln(a/a*o)] + sinh^ 77*) ^a^/m. (79) 

The coefficient of ha^/m is shown as a function of a/a* in Fig. 7 for several 
values of 77*. As a varies, the coefficient of ha^/m oscillates between about 
67.12 e"^''* sinh^ 77* and about 67.12 e^^'^* cosh^ 7/*. Thus one effect of the deep 
dimers is to eliminate the zeros in ccshaiiow at a = (e'^/*°)"a*o- The depth of 
the minimum is quadratic in 77* as 77* — > 0, so the coefficient of ha^/m at 
a — (e'^/*°)"'a*o can be very small if the Efimov width parameter 77* is small. 

The existence of deep dimers opens up additional channels for 3-body recom- 
bination. If there are no deep dimers, 3-body recombination can only produce 
the shallow dimer if a > and it cannot proceed at all if a < 0. If there are 
deep dimers, they can be produced by 3-body recombination for either sign 
of a. We denote by ctdeep the inclusive contribution to the event rate constant 
defined in Eq. (62) from 3-body recombination into all the deep dimers. 

If a > 0, the analytic expression for ctdeep is 



^ _ 647r2(47r-3V3)coth(7rso)sinh(277*) ^a^ ^^^^ 
'^^^'^ sinh^(7rso + ?7*) + cos2[soln(a/a*o)] "t. 

where a*o ~ 4.5 a* is given in Eq. (64). The coefficient of Ua^/m has very weak 
log-periodic dependence on are*. We can use the appproximation sinh(7rso + 
77*) fa e^* sinh(7rso) to simphfy the expression in Eq. (80): 



'-'^deep 



16.84 (1 - e-^"*) ha^/m , (a > 0) . (81) 



The numerical result for the coefficient in Eq. (81) was first derived in Ref. [38]. 
The coefficient of ha'^/m, which is independent of a/a*, is shown in Fig. 7 for 
several values of 77*. In the limit 77* — > 0, ctdeep approaches zero linearly in 77*. 
In the limit 77* — > 00, the rates a^eep in Eq. (80) and ccghaiiow in Eq. (78) are 
almost equal, differing only by the factor coth(7rso) ~ 1.004. 

If a < 0, the 3-body recombination rate constant is 



^ _ 4590sinh(277*) ha"^ (a < 0) (82) 

'^^^^ sin^[so ln(a/a'^)] + sinh^ 7^* m ' ' 
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Fig. 8. The 3-body recombination rate constant Odeep in units of ha^/m as a function 
of a/ a'^ for a < and several values of . 

where a'^ ^ — 1.6^"^ is given in Eq. (56). The explicit formula for adeep in 
Eq. (82) was first derived in Ref. [38] [£j The coefficient 4590 is only known 
numerically. The coefficient of ha^/m is shown as a function of a/a'^ in Fig. 8 
for several values of rj^. It displays resonant behavior with maxima when the 
scattering length has one of the values (e'^/*°)"a'^ for which there is an Efimov 
trimer near the 3- atom threshold. In the limit 77* — 0, the maximum value 
9180coth?7* of the coefficient of %a^/m in Eq. (82) diverges. In the limit 77* — > 
00, the log-periodic dependence of the coefficient on an^ disappears and it 
approaches the constant 9180. 



6 Applications to Atoms 



In this section, we describe applications of universality to '^He atoms and to 
alkali atoms near a Feshbach resonance. 



6.1 Helium atoms 



Helium atoms provide a beautiful illustration of universality in the 3-body 
system [36]. The interatomic potential between two ^He atoms does not sup- 
port any deep dimers, so universality is realized in its simplest form with a 
single 3-body parameter k^,. The binding energies of the ^He trimers have 
been calculated accurately for a number of different model potentials for the 
interaction between two ^He atoms. For the purposes of illustration, we will 

^ See also Ref. [40] for an earlier calculation that was perturbative in ry*. 
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use the TTY potential [9]. The scattering length for the TTY potential is 
a = 188.99 flQ. This is much larger than its effective range = 13.85 a^, which 
is comparable to the van der Waals length scale £vdw = 10.2 oq. The TTY 
potential supports a single 2-body bound state, the ^He dimer whose binding 
energy is E2 — 1.30962 mK. This binding energy is small compared to the 
natural low-energy scale for ^He atoms: Ey^w ~ 420 mK. The TTY potential 
supports exactly two 3-body bound states: the ground-state trimer, which we 
label n = 0, and the excited trimer, which we label n = 1. There have been 
several accurate calculations of the binding energies eP and E^^^ for the TTY 
potential [41,42,43]. The results agree to within 0.5% for both E^^^ and Ei^\ 
The results of Ref. [42] are ^ = 126.4 mK and E^^^ = 2.28 mK. 

The ^He dimer was first observed in 1992 by the Minnesota group [44]. They 
used the expansion of ''He gas at room temperature from a pulsed valve into 
a vacuum chamber to create a beam of ^He atoms and clusters with a trans- 
lational temperature near 1 mK. The dimers were detected by using electron 
impact ionization to produce helium dimer ions that were observed by mass 
spectrometry. In another experiment in 1995, they determined the size of the 
^He dimer by measuring the relative transmission rates of He atoms and dimers 
through a set of nanoscale sieves [45] . Their result for the mean separation of 
the atoms in the dimer was (r) = (62 ± 10) A. 

The ''He dimer was also detected nondestructively by the Gottingen group 
[46] . They made a beam of ^He atoms and clusters with temperature ranging 
from 6 K to 60 K by allowing cryogenic ^He gas to escape from an orifice. The 
beam was passed through a nanoscale transmission grating, and the dimers 
were detected by observing a diffraction peak at the expected angle. In 1995, 
this experiment was used to make the first observation of the ground state "^He 
trimer [4] . The excited ^He trimer has not yet been observed. In a subsequent 
experiment, the ^He tetramer was observed and the formation rates of dimers, 
trimers, and tretramers was measured as functions of the temperature and 
pressure of the gas from which the beam escaped [47]. In 2000, the Gdttingen 
group determined the size of the dimer from the relative strengths of the 
higher order diffraction peaks from the nanoscale transmission grating [48]. 
Their result for the mean separation of the atoms is (r) = (52 ± 4) A. This 
result is in good agreement with the theoretical prediction from the TTY 
potential: (r) = a/2 = 50 A. From this measurement, the binding energy of 
the ^He dimer was inferred to be E2 = (l.ltjj ^) The scattering length for 
He atoms was inferred to be a = (104t^8) A. 

Lim, Duffy, and Damert proposed in 1977 that the excited state of the ^He 
trimer is an Efimov state [49]. This interpretation is almost universally ac- 
cepted. Some researchers have proposed that the ground state trimer is also 
an Efimov state [18,50,51]. This raises an obvious question: what is the def- 
inition of an Efimov state? The most commonly used definition is based on 
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rescaling the depth of the 2-body potential: V{r) — > AV^(r). According to 
the traditional definition, a trimer is an Efimov state if its binding energy as 
a function of the scaling parameter A has the qualitative behavior illustrated 
in Fig. 4. As A is decreased below 1, the trimer eventually disappears through 
the 3-atom threshold. As A is increased above 1, the trimer eventually disap- 
pears through the atom-dimer threshold. Calculations of the trimer binding 
energies [52] using a modern helium potential show that the excited trimer 
satisfies this definition of an Efimov state but the ground state trimer does 
not. The excited trimer disappears through the 3-atom threshold when A is 
decreased to about 0.97, and it disappears through the atom-dimer threshold 
when A is increased to about 1.1. The ground state trimer disappears through 
the 3-atom threshold when A is about 0.9. However, as A is increased above 1, 
its binding energy relative to the atom-dimer threshold continues to increase. 
Thus the ground state ^He trimer does not qualify as an Efimov state by the 
traditional definition. 

The traditional definition of an Efimov state described above is not natural 
from the point of view of universality. The essence of universality concerns the 
behavior of a system when the scattering length becomes increasingly large. 
The focus of the traditional definition is on the endpoints of the binding energy 
curve in Fig. 4, which concerns the behavior of the system as the scattering 
length decreases in magnitude. The problem is that the rescaling of the poten- 
tial can move the system outside the universality region defined by \a\ ^ 
before the trimer reaches the atom-dimer threshold. We therefore propose a 
definition of an Efimov state that is more natural from the universality per- 
spective. A trimer is defined to be an Efimov state if a deformation that tunes 
the scattering length to ±oo moves its binding energy along the universal 
curve illustrated in Fig. 4. The focus of this definition is on the resonant limit 
where the binding energy crosses the 1/a = axis. In particular, the binding 
energy at this point should be larger than that of the next shallowest trimer 
by about a factor of 515. In the case of ^He atoms, the resonant limit can 
be reached by rescaling the 2-body potential by a factor A ~ 0.97 [52]. At 
this point, the binding energy of the ground state trimer is larger than that 
of the excited trimer by about a factor of 570. The closeness of this ratio to 
the asymptotic value 515 supports the hypothesis that the ground state ^He 
trimer is an Efimov state. 

In order to apply the universal predictions for low-energy 3-body observables 
to the case of "^He atoms, we need a 2-body input and a 3-body input to 
determine the parameters a and k*. The scattering length a — 188.99 for the 
TTY potential can be taken as the 2-body input. An alternative 2-body input 
is the dimcr binding energy: E2 = 1.31 mK. A scattering length can be 
determined by identifying E2 with the universal binding energy of the shallow 
dimer: E2 = h?/{ma%). The result for the TTY potential is an = 181.79 ao. 
The 3.8% difference between a and a£) is a measure of how close the system is 
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to the scaling limit. To minimize errors associated with the deviations of the 
system from the scaling limit, it is best to take the shallowest 3-body binding 
energy available as the input for determining In the case of *^He atoms, 
this is the binding energy E^^^ of the excited trimer. Experience has shown 
that the universal predictions are considerably more accurate if E2 and £'3^^ 
are used as the inputs instead of a and £'3^^ [36] . 

We proceed to consider the universal predictions for the trimer binding ener- 
gies. Having identified E^^ with the universal trimer binding energy , we 
can use Efimov's binding energy equation (44) with = 1 to calculate up 
to multiplicative factors of e'^/'o ^ 22.7 [36]. The resuh is = 0.00215 or 
= 0.00232 Cq^, depending on whether E2 or a is used as the 2-body input. 
The intuitive interpretation of k^, is that if a parameter in the short-distance 
potential is adjusted to tune a to -|-oo, the binding energy Ei^^ should ap- 
proach a limiting value of approximately h^nl/m, which is 0.201 mK or 0.233 
mK depending on whether the 2-body input is E2 or a. 
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TTY potential 


100.0 


96.2 


2.28 


126.4 




universality 




input 


input 


129.1 


5.38 X 10^ 


universality 


input 




input 


146.4 


6.23 X 10^ 



Table 1 

Binding energies of the ^He trimers for the TTY potential (row 1) compared 
to the universality predictions using as the inputs cither E2 and E^^ (row 2) or a 
and £'3^'* (row 3). Lengths are given in A and energies are given in mK. The trimer 
binding energies for the TTY potential are from Ref. [42]. (Note that h'^/m = 
12.1194 for ^He atoms.) 

Once has been calculated, we can solve Eq. (44) for the binding energies 
of the deeper Efimov states. The prediction for the next two binding energies 
are shown in Table 1. The prediction for E^^ differs from the binding energy 
of the ground-state trimer by 2.6% or 16.4%, depending on whether E2 or a 
is taken as the 2-body input. The expected order of magnitude of the error is 
the larger of -^vdw/o = 5.4% and {E^f ^/EydwY^^ ~ 50%. The errors are much 
smaller than (4'^^/E'vdw)^^^, suggesting that the scahng hmit is more robust 
than one might naively expect. Efimov's equation (44) also predicts infinitely 
many deeper 3-body bound states. The predictions for the binding energy 
Ejf^^ of the next deepest state are given in Table 1. The predictions are more 
than two orders of magnitude larger than the van der Waals E'vdw ~ 420 mK. 
We conclude that this state and all the deeper bound states are artifacts of 
the scaling limit. 

Using the above values of k,^ for the TTY potential, we can immediately predict 
the atom-dimer scattering length a ad- If a is used as the 2-body input, we 
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find a AD ~ 0.94 a, corresponding to a ad ~ 178 ao- If E2 is used as the 2-body 
input, we find a ad ~ 1.19 a^, corresponding to a ad ~ 216 Oq. These values are 
in reasonable agreement with the calculation of Ref. [43], which gave a ad ~ 
248(10) Oq. Since r^/a = 7.3% for the TTY potential, much of the remaining 
discrepancy can perhaps be attributed to effective-range corrections. 

Universality can also be used to predict the 3-body recombination rate con- 
stant for *^He atoms interacting through the TTY potential. The prediction for 
Qishaiiow is 2.9 ha% /m or 6.9 Ua^ / m, depending on whether E2 or a is used as the 
2-body input. In either case, the coefficient of ha^/m is much smaller than the 
maximum possible value 67.1. Thus ^He atoms are fortuitously close to a com- 
bination of a and /t* for which ashaiiow is zero. The 3-body recombination rate 
constant has not yet been calculated for the TTY potential. It has been calcu- 
lated for the HFD-B3-FCI1 potential [53]. The result is ctshaUow = 12 x lO^^*^ 
cm^/s [37], which compares well with the universal prediction 9 x 10~^^ cm^/s 
obtained using E2 and £'3^-' as the inputs [6] . 



6.2 Alkali atoms near a Feshbach resonance 

Alkali atoms near a Feshbach resonance provide a unique window on Efimov 
physics, because the scattering length can be tuned experimentally. Inelastic 
loss rates have proved to be an especially powerful probe of 3-body processes 
in these systems. In this subsection, we discuss some key experiments on 3- 
body losses for ^^Na, ^^Rb, ^^Rb, and ^^^Cs atoms near Feshbach resonances 
and compare them with the predictions from universality. 

The first observation of the enhancement of inelastic losses near a Feshbach 
resonance was by the MIT group [14]. They created Bose- Einstein condensate 
of ^^Na atoms in the |1, -|-1) hyperfine state and used a magnetic field to ad- 
just the scattering length. They observed enhanced losses near the Feshbach 
resonances at 907 G and 853 G. Since the |1, +1) state is the lowest hyperfine 
state, the inelastic losses come primarily from 3-body recombination. The 3- 
body losses near the Feshbach resonance at 907 G were studied systematically 
by the MIT group [15]. They also studied 3-body losses in a Bose-Einstein 
condensate of ^'^Na atoms in the |1, —1) hyperfine state near a Feshbach res- 
onance at 1195 G. The off-resonant scattering length in this region of high 
magnetic field is approximately 52 ao, which is significantly smaller than the 
van der Waals scale ^vdw ~ 90 oq. Near the Feshbach resonances at 907 G, 
the authors of Ref. [15] were able to increase the scattering length by about 
a factor of 5 over the off-resonant value. However this it still not very large 
compared to the i^dw, so the universal theory does not apply. A quantitative 
description of the data on 3-body recombination near the 907 G resonance 
has been given using a scattering model for a Feshbach resonance with zero 
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off-resonant scattering length [54]. 




Fig. 9. The 3-body loss coefficient K3 in Rb as a function of the magnetic field 
near the Feshbach resonance at B = 155 G. The dots are the measured values of 
from Ref. [55], while the triangles indicate upper bounds on K^. The curves are 
for K*ao = 2 and several values of 

The inelastic collision rate for ultracold ^^Rb atoms in the |2, —2) hyperfine 
state near the Feshbach resonance at 155 G has been studied by the JILA 
group [55]. By exploiting the different dependences on the number density, 
they were able to separate the contributions from 2-body and 3-body pro- 
cesses. They measured the 2-body and 3-body loss coefficients as a function 
of the magnetic field B from 110 G to 150 G. At some values of B, they were 
only able to obtain an upper bound on the 3-body loss coefficient K3 = 3a. 
Their results for near the Feshbach resonance at 155 G are shown in Fig. 9. 
The universal results for the 3-body recombination rate at threshold are given 
by Eq. (82) for a < and they can be approximated by the sum of Eqs. (79) 
and (81) for a > 0. A precise determination of the Efimov parameters and 
?7* is not possible with the data in Fig. 9. However as shown in Fig. 9, the data 
can be described reasonably well by the universal formulas with the values 
K^,ao = 2 and rj^ = 0.02. The curves for different values of rj^ illustrate how the 
minima in the recombination rate for a > are filled by recombination into 
deep bound states. 

The 3-body losses in ultracold ®^Rb were studied by the Garching group [56]. 
Roughly 90% of the atoms were in the |1,+1) hyperfine state. Most of the 
remaining atoms were in the |1,0) state, but there were also some atoms in 
the 1 1,-1) state. By monitoring the atom loss, the authors observed more 
than 40 Feshbach resonances for various combinations of hyperfine states of 
^^Rb at magnetic fields ranging from 300 G to 1300 G. Measurements of the 
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resonances were used to deduce an improved atomic potential for ^^Rb. Away 
from the Fcshbach resonances, the scattering lengths for ^''Rb have natural 
values comparable to the van der Waals scale ^vdw ~ 165 Oq. The 3-body loss 
rate — Sa in the vicinity of the resonance at 1007 G was measured as a 
function of the magnetic field. Away from the resonance, K3 was measured to 
be 3.2(1.6) X 10~^^ cm^/s. Near the resonance, the rate constant was observed 
to increase by as much as a factor of 300. 

The 3-body recombination rate in a Bose-Einstein condensate of ^^Rb atoms 
in the |1,+1) hyperfine state has been measured just below the Feshbach 
resonance at 1007 G by the Garching group and by the Oxford group [57]. The 
values of the magnetic fields ranged from about 1 G below the resonance, where 
the scattering length has a natural value comparable to the van der Waals scale 
^vdw ~ 165 Oo, to about 0.03 G below the resonance, where the scattering 
length is about a factor 5 larger than £vdw- Only the last few data points are 
in the universal region. The measured rate constant = 3a appears to scale 
as a smaller power of a than the scaling prediction a'^. The measurements of 
agree reasonably well with the results of exact solutions to the 3-body 
Schrodinger equation for a coupled-channel model [57]. These calculations 
predict a local minimum of K3 that can be attributed to Efimov physics at 
a magnetic field about 0.015 G below the resonance, just outside the range 
of the experiments. The atom-dimer relaxation rate for the coupled-channel 
model was also calculated in Ref. [57]. The model predicts a resonance in the 
relaxation rate less than 0.1 G below the resonance. 

The 3-body recombination rate for an ultracold gas of ^^^Cs atoms in the 
|3, +3) hyperfine state was measured as a function of the magnetic field by 
the Innsbruck group [58]. The 3-body recombination loss rate = 3a is 
very small near 17 G, where the scattering length goes through a zero that 
comes from the interplay between a broad Feshbach resonance and the large 
off-resonant scattering length for ^^^Cs. For magnetic fields between 50 G 
and 150 G, the measurements confirmed the power-law dependence of 
predicted by scaling. In this region of the magnetic field, the scattering length 
increased from about 1000 Qq to about 1600 Gq. Since the scattering length 
increased only by a factor of 1.6, it was not possible to observe any logarithmic 
variations of the coefficient of a^. 

The Innsbruck group extended their measurement of the 3-body recombination 
loss rate = 3a for ^'^^Cs atoms in the |3, 4-3) hyperfine state to values of 
the magnetic field below 17 G where the scattering length is negative [5]. The 
negative scattering length could be made as large as — 2500ao. They observed 
a resonant enhancement of for a near — SOOao- The enhancement can be 
explained by the presence of an Efimov state near the 3-atom threshold. At the 
lowest temperature they were able to reach, which was 10 nK, the shape of the 
resonance could be fit reasonably well using the universal result in Eq. (82), 
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Fig. 10. The 3-body loss coefficient in "'^Cs for negative values of the scattering 
length near the Efimov resonance at a ~ — 850ao. The data points are for T = 10 
nK and are taken from Ref. [5]. The curves are for k^, = 0.945/ao and three different 
values of r/* . 

which has a peak near a*. The best fit values of the Efimov parameters are 
a'^ = —850(20) Qq and r/* = 0.06(1). In Fig. 10, we compare the data of Ref. [5] 
with the universal result. The parameters a'^ = — 850ao (which corresponds 
to = 0.945 flg^) and ?7* = 0.06 give a good fit to the data. Both Efimov 
parameters are well determined by the data: k^, by the position of the Efimov 
resonance and 77* by the height and width of the peak. Since the van der Waals 
length scale for ^^^Cs is £vdw ~ 200 Oq, the resonance is reasonably deep into 
the universal region. Thus these results from Innsbruck group seem to be the 
first experimental evidence for the existence of Efimov states [5] . 

The Innsbruck group also carried out a more careful analyisis of 3-body re- 
combination at magnetic fields above 170 G, where a is positive [5]. They 
measured the loss rate with higher accuracy and at more values of the 
magnetic field. They observed a local minimum in the atom loss rate for a 
near 210 Oq. The dependence of the loss rate constant K3 on the scattering 
length can be fit reasonably well with the universal formula given by the sum 
of Eqs. (79) and (81). The value of 77* is not well determined by the data and 
the fit yields only the upper bound 77=,, < 0.2. The best fit for the value of a at 
which the coefficient of is maximal is a+ = 1060(70) a^. This corresponds 
to a* ~ 1170 Oq. The positive scattering length data from Ref. [5] are shown in 
Fig. 11. The Efimov parameter = 0.707ag^ corresponding to = 1170 Oq 
gives a good fit to the data above a ~ 400 oq. If the scattering length could 
be increased past cxd to a region where a is large and negative, universality 
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Fig. 11. The 3-body loss coefficient in ^^■^Cs for positive values of the scattering 
length. The data points are for T = 200 nK (dots) and for temperatures between 
250 nK and 400 nK (triangles) and they are taken from Ref. [5] . The curves are for 
= 0.707/ao and three different values of r|^,. 

predicts that there should be an Efimov state at the 3-atom threshold when a 
is approximately —1100 qq. This is not far (at least on a log scale) from the po- 
sition —850 Oq where the resonance was observed. However the region of large 
negative a in which the resonance is observed is separated from the region of 
large positive a by a region where a passes through zero and universal theory 
does not apply. As a consequence, the Efimov parameters for the two regions 
need not be equal. One note of caution in applying the universal expression 
for the 3-body recombination rate at threshold is that the positive scattering 
length data in Fig. 11 was obtained at temperatures that ranged from 200 nK 
to 400 nK, which is much larger than the temperature 10 nK for the negative 
scattering length data in Fig. 10. The temperature for the positive scattering 
length data may not be low enough to use the universal expression for the 
3-body recombination rate at threshold. 

The Innsbruck group has also observed inelastic losses of shallow dimers com- 
posed of ^^^Cs atoms in the |3, +3) hyperfine state [59]. The fraction of dimers 
that are lost remains roughly constant for magnetic fields in the range from 
between 14 G and 19.8 G. However, the fraction increases rapidly toward 1 
as the magnetic field increases further toward a Feshbach resonance at 19.84 
G. In this experiment, resonant enhancements of the inelastic loss rate were 
also observed at lower values of the magnetic field. They can be attributed 
to effects of a ^^^Cs tetramer near the threshold for two ^^^Cs dimers. In this 
region of the magnetic field, the dimers are deeply bound, so these results have 
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nothing to do with universahty. 



Very recently, the Innsbruck group has created mixtures of ^^^Cs atoms in the 
1 3, +3) hyperfine state and shallow dimers composed of those atoms [60]. They 
have observed a resonance enhancement in the inelastic atom-dimer collision 
rate near a ^ 400 oq [60]. The enhancement can be explained by the presence 
of an Efimov state near the atom-dimer threshold. They have measured the 
rate constant (3 defined by Eq. (75) as a function of the magnetic field at 
a temperature near 250 nK. The universal predictions for /3 as a function 
of T have been calculated for temperatures small compared to the binding 
energy of the shallow dimcr. It might be possible to use these measurements 
to determine the Efimov parameters a* and 7]^ accurately for this region of 
large positive scattering length [61]. 



7 Efimov Physics in Other Three-body Systems 

Very few universal results have been calculated for 3-body systems other than 
identical bosons. They are summarized in Ref. [6]. In this section, we consider 
only the basic question of whether the Efimov effect occurs in the 3-body 
system. This question can be answered by determining the channel eigenvalue 
Ao(i?) for the lowest hyperspherical potential in the scaling limit. The Efimov 
effect occurs if \o{R) is negative at R = 0. If Ao(0) = — Sq, the discrete scaling 
factor for Efimov physics is e^^^°. 

7.1 Unequal scattering lengths 

There are 3-body systems in atomic physics in which the three atoms all have 
the same mass, but the three pairs of atoms need not all have the same scat- 
tering lengths. For example, different hyperfine spin states of the same atom 
have the same mass, but the scattering lengths be different for each 

pair ij of hyperfine states. As another example, different isotopes of a heavy 
atom have nearly the same masses. It is therefore worthwhile to consider the 
universal behavior of systems with equal masses and with scattering lengths 
that are large but not necessarily equal. 

The Efimov effect in general 3-body systems was first discussed by Amado and 
Noble [62] and by Efimov [63,64]. A summary of their results is as follows. If 
only one of the three scattering lengths is large, the Efimov effect does not 
occur. If two of the scattering lengths are large, the Efimov effect occurs with 
a discrete scaling factor of 1986.1 unless two of the three particles are identical 
fermions, in which case the Efimov effect does not occur. If all three scattering 
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lengths arc large, the Efimov effect occurs with a discrete scaling factor of 22.7. 
For a derivation of these results, see Ref. [6]. 

If the atoms are distinguishable but related by an internal symmetry, the 
condition for the Efimov effect is modified. An example is the spin states 
belonging to a hyperfine multiplct at zero magnetic field, which arc related by 
a spin symmetry. A general treatment of the effects of an internal symmetry 
was given by Bulgac and Efimov in Ref. [65]. 



1.2 Unequal masses 



In the 2-body sector, the universal results for particles of unequal masses 
are only a little more complicated than those for the equal-mass case in Sec- 
tions 3.1 and 3.2. Let the atoms 1 and 2 have a large scattering length 
and unequal masses mi and m2. If ai2 is large and positive, the atoms 1 and 
2 form a shallow dimer with binding energy 
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where mi2 = mim2/(mi -|- 7712) is the reduced mass. 



In the general 3-body system, the three masses can be unequal and any com- 
bination of the three scattering lengths can be large. The Efimov effect in 
general 3-body systems was first discussed by Amado and Noble [62] and by 
Efimov [63,64]. The special case in which two of the three particles have the 
same mass was also discussed by Ovchinnikov and Sigal [66]. The conditions 
for the existence of the Efimov effect and the value of the discrete scaling 
factor depend on the ratios of the masses. We first summarize the results for 
the extreme cases in which two masses are equal and the third mass is either 
much larger or much smaller. In the case of two heavy particles and one light 
particle, the Efimov effect occurs provided the heavy-light scattering length is 
large. The Efimov effect can be understood intuitively in this case by using the 
Born-Oppenheimer approximation [67] . In the case of one heavy particle and 
two light particles, the Efimov effect occurs only if all three pairs of particles 
have large scattering lengths. 

We next consider a more general case in which all three pairs of atoms have 
a large scattering length. This excludes the possibility of any pair of particles 
being identical fermions. The condition for the Efimov effect does not depend 
on the values of the scattering lengths as long as their absolute values are large 
compared to the range. The Efimov effect occurs for any values of the masses. 
The discrete scahng factor is largest if all three masses are equal. It has the 
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Fig. 12. Discrete scaling factor e^l^^ for two particles of equal mass m\ = m2 
as a function of the mass ratio nii/m^ for the case in which all three pairs have 
large scattering lengths. The particles 1 and 2 can be either identical bosons or 
distinguishable. The dot indicates the case of three identical bosons. 

same value e'^/*" ~ 22.7 as for three identical bosons. In the case of two equal- 
mass particles, the discrete scaling factor is the same whether the equal-mass 
particles are identical bosons or distinguishable. The discrete scaling factor for 
mi = m2 is shown as a function of the mass ratio mi/m^ = m2/m^ in Fig. 12. 
In the limit mi = m2 ^ m^, the discrete scaling factor approaches 15.7. In 
the limit mi = m2 ^ m^, it approaches 1. 
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Fig. 13. Discrete scaling factor e^^^'^ for two particles of equal mass rrii = m2 as a 
function of the mass ratio nii/m^ for the cases in which two pairs have large scatter- 
ing lengths. If 023 and 031 are large, particles 1 and 2 can be either identical bosons 
or distinguishable particles (solid line) or else identical fermions (dash-dotted line). 
If ai2 and 031 (or ai2 and 023) are large, particles 1 and 2 must be distinguishable 
particles (dashed line). 

We now consider the more general case in which only two pairs of atoms have 
large scattering lengths. For simplicity, we consider only the special cases in 
which particles 1 and 2 have the same masses mi = m2. There are three cases 
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to consider. The discrete scaling factors are shown in Fig. 13 as a function of 
the mass ratio mi/m^ for all three cases. If the large scattering lengths are 023 
and aai and if particles 1 and 2 are either identical bosons or distinguishable, 
the discrete scahng factor e'^/*° decreases monotonically from 00 to 1986. If the 
large scattering lengths are 023 and 031 and if particles 1 and 2 are identical 
fermions, the Efimov effect occurs only if mi/ms exceeds the critical value 
13.607. The discrete scaling factor decreases monotonically from 00 to 1 as 
TTii/ms increases from the critical value to 00. If the large scattering lengths 
are ai2 and 031 (or ai2 and 023), particles 1 and 2 must be distinguishable. 
The discrete scaling factor increases monotonically from 94.36 to 1986.1 to 00 
as mi/ms increases from to 1 to 00. 



8 Extensions and Outlook 

In this section, we discuss extensions of universality beyond the 3-body prob- 
lem in the scaling limit: effective-range corrections, microscopic models with 
a large scattering length, and the 4-body problem. We end this review with 
an outlook on the study of Efimov physics. 



8.1 Effective range corrections 

The scaling limit of a system with a large scattering length is the limit in which 
the range of the interaction is taken to zero. There are corrections to the scaling 
limit from the fact that the range for any physical system is nonzero. These 
corrections can be calculated as an expansion in powers of £/a, where i is the 
natural low-energy length scale. The most important corrections come from 
the effective range r^, which is defined by the expansion of the S-wave phase 
shift in Eq. (5). Examples of effective range corrections for 2-body observables 
are given in Eqs. (33) and (34). 

The most important correction to the scaling limit in the 3-body system also 
comes from the effective range r^. At first order in the expansion in £/\a\, 
the effective range is the only new parameter that enters. Efimov has shown 
that the effective range leads to a correction term r^/ in the scale-invariant 
hyperspherical potential in Eq. (21) [68]. He gave expressions for the first-order 
effective range corrections to the binding energy of the Efimov trimer and to 
the atom-dimer scattering length [69]. Moreover, he applied his results to the 
3-nucleon system, calculating the range corrections to the correlation between 
the triton binding energy and the spin-doublet neutron-deuteron scattering 
length [70] and to the spin-quartet neutron-deuteron scattering length [69] . 
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Here, we will give only a brief summary of the systematics of higher-order range 
corrections in channels where the Efimov effect occurs. For a more complete 
discussion, we refer the reader to our longer review [6]. Since the correction 
to the universal long-range potential is singular, the range corrections require 
renormalization. Effective field theory provides a powerful framework to deal 
with these naively divergent corrections. The first-order effective range cor- 
rection to the spin-doublet neutron-deuteron scattering phase shifts was first 
calculated using this method in [71] and later confirmed using a different renor- 
malization scheme [72] . The systematics of contributions beyond first order in 
l/\a\ was studied by various authors [73,74,75,76,77]. The first calculation of 
the triton channel of the 3-nucleon system at second order in (■/\a\ was carried 
out in Ref. [73]. As suggested by naive dimensional analysis, a new 3-body 
parameter was included in this calculation at order (£/|a|)^. Recently, Plat- 
ter and Phillips revisited this problem using a subtractive renormalization 
scheme. They found that the scattering amplitudes could be renormalized at 
order {C./\a\Y without the need of a new subtraction constant [76,77]. Accord- 
ing to their analysis, a new 3-body parameter is not required until at least 
order (£/|a|)^. They applied their method to the system of "^He atoms and 
found good convergence of the expansion in £/|a| [76]. 



8.2 Microsopic models 

Since the low-energy behavior of particles with large scattering lengths is in- 
sensitive to the details of their interactions at short distances, the potential 
provides an inefficient encoding of the relevant physics. The sensitivity to short 
distances enters primarily through the scattering length and other constants 
that describe low-energy scattering. This motivates a more phenomenological 
approach in which the interactions are described by a scattering model, which 
can be specified by a parameterization of low-energy scattering amplitudes. 
The parameters of the scattering model can be treated as phenomenological 
parameters that can be tuned to reproduce the observed low-energy observ- 
ables of the 2-body system. Scattering models will exhibit universal low-energy 
behavior as the scattering length a is tuned to ±oo. Low-energy 3-body ob- 
servables reduce in this limit to functions of a and the Efimov parameters 
and T]^, all of which are functions of the parameters of the scattering model. 

An example of a scattering model that may provide an accurate description 
of atoms near a Feshbach resonance is the resonance model [78]. This model 

provides a natural description of atoms that have a weak coupling to a closed 
channel in which there is a diatomic molecule extremely close to the atom- 
atom threshold. The resonance model has three parameters that can be defined 
by specifying the S-wave phase shift 5o{k): 
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kcotdo{k) = — A + 




(84) 



The scattering length and the effective range are a = X — g'^/v and = 
—^{g — i/X/gY, respectively. Thus the effective range is negative definite. The 
approximation in Eq. (8) for the scattering length as a function of the magnetic 
field B for an atom near a Feshbach resonance can be reproduced by taking 
A and g to be independent of B and v to be linear in B — i?rcs- Thus A can 
be identified with the off-resonant scattering length, v is proportional to the 
detuning energy of a molecule in a closed channel, and g is the strength of 
the coupling to the closed channel. The resonance model will exhibit universal 
behavior as v is tuned to zero. Three-body obscrvables will approach universal 
functions of a, k*, and 7]^. The Efimov parameter must be a function of A 
and g, while 7]^ must be a function of the dimensionless combination \g^. 

Some 3-body observables have been calculated in the resonance model in the 
special case A = [54]. The scattering length and the effective range are 
a = —g^/i^ and Vg = —2/g^, respectively. It is convenient to define a param- 
eter i?* = 1/g'^, with dimension of length. The atom-dimer scattering length 
qad and the 3-body recombination rate into the shallow dimer a have been 
calculated in this model as functions of a and R^, [54]. For a ^ R^, they have 
the universal behavior given in Eqs. (61) and (63). The 3-body parameter 
was determined in Ref. [54] to be 

So ln(«;^) = So ln(2.5/i?*) mod tt . (85) 

The cross-over to the universal behavior occurs when a is comparable to i?*. 
The first divergence of the atom-dimer scattering length occurs at a = 0.45 R^. 
This marks the emergence of the first Efimov state below the atom-dimer 
threshold. The second and higher divergences of gad occur at values of a 
that are well-approximated by the universal predictions: a — (e'^/**'')" 0.64 i?*, 
n — 1,2,..., where e^^^° ~ 22.7. The first zero of the recombination rate 
constant occurs at a = 3.3 i?*. The second and higher zeroes occur at values of 
a that are well approximated by the universal predictions: a — (e^'/^o)" 2.9 i?*. 



8.3 Four-body problem 

In the 3-body problem, exact numerical solutions are facilitated by the Fad- 

dccv equations. The generalization of the Faddeev equations to the A^-body 
problem with N > 4 was given by Yakubovsky [79] . An equivalent set of equa- 
tions was given independently by Grassberger and Sandhas [80]. Due to the 
complexity of these equations, exact numerical solutions ior N — A have only 



n = l,2,.... 
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recently been obtained. In atomic physics, the only numerically exact A^-atom 
calculations for > 4 that we are aware of arc for ground-state binding ener- 
gies. There are no such calculations for weakly-bound A^-atom molecules that 
might be governed by universality. 

Testing universality in systems with four and more particles is a new frontier. 
Even theoretically, the universal properties in the 4-body sector arc not fully 
understood. An important issue in the 4-body system with a large scattering 
length is how many parameters are required to describe the system in the scal- 
ing limit, that is, up to corrections that decrease like ^/\a\ asa— >±oo. In the 
case of identical bosons, low-energy 4-body observables necessarily depend on 
the 2-body parameter a and the 3-body parameter k^,. But are any additional 
4-body parameters required? There are theoretical arguments in support of 
both answers to this question. There is a renormalization argument for zero- 
range 2-body potentials that indicates that an additional A'^-body parameter is 
required to calculate A'^-body binding energies for all A^ > 3 [81]. On the other 
hand, a power-counting argument within the effective field theory framework 
suggests that after adding the 3-body parameter k*, no additional parameters 
are necessary to calculate A'"-body observables for A?" > 4 [82] . 

Platter et al. have recently studied the universal properties of the system of 
four identical bosons with short-range interactions in an effective quantum 
mechanics approach [83]. They constructed the effective interaction poten- 
tial at leading order in the large scattering length and computed the 4-body 
binding energies using the Yakubovsky equations. They found that cutoff in- 
dependence of the 4-body binding energies does not require the introduction 
of a 4-body parameter at leading order in l/\a\. This suggests that 2-body 
and 3-body interactions arc sufficient to rcnormalizc the 4-body system. They 
applied their equations to '^He atoms and calculated the binding energies of 
the ground state and the excited state of the ^He tetramer. Using the bind- 
ing energy E2 of the ^He dimer as the 2-body input and the binding energy 
£^3 of the excited state of the '^Hc trimcr as the 3-body input, they found 
good agreement with results of Blume and Greene [84] calculated using an 
approximate numerical method that combines Monte Carlo methods with the 
hyperspherical expansion. 

Yamashita et al. [85] recently claimed that for a 4-body system close to a 
Feshbach resonance, a new 4-body parameter will enter already at leading 
order in ^/\a\. They motivated their claim from a model-space reduction of a 
realistic 2-body interaction. In the framework of the renormalized zero-range 
model, they found a strong sensitivity of the 4-body ground state energy 
to a 4-body subtraction constant in their equations. The results of Ref. [83] 
for the special case of "^He atoms were also reproduced. A drawback of their 
analysis is the focus on the deepest 4-body state only. It remains to be seen 
whether their findings are universal or just artefacts of their regularization of 
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the renormalized zero-range model. 



8.4 Outlook 

Vitaly Efimov's discovery in 1970 of the Efimov effect in identical bosons was 
like finding a nugget of a precious metal. Other physicists quickly verified the 
amazing properties of this nugget and tried to find similar nuggets in other 
systems. Efimov went on to show that the nugget came from a rich vein of ore 
corresponding to universal properties of 3-body systems with a large scattering 
length. He pointed out that the Efimov effect was just a hint of a discrete 
scaling symmetry, and showed how the discrete scaling symmetry could be 
exploited to mine the ore. However not many physical systems in which Efimov 
physics plays an important role were found. In those that were identified, such 
as ^He atoms, it was difficult to study Efimov physics experimentally. Efimov's 
mine was eventually closed due to lack of interest from other physicists as well 
as the lack of a market for the ore. 

The new field of cold atom physics has changed the situation dramatically. The 
development of the technology for trapping alkali atoms and cooling them to 
ultralow temperatures has provided a long list of 3-body systems in which 
Efimov physics can be studied. Especially important is the possibility of using 
Feshbach resonances to control the scattering length and make it arbitrarily 
large. This development has justified the reopening of Efimov's mine. 

The theoretical study of Efimov physics is still in its infancy. In the case of 
3 identical bosons, the rates for most of the basic scattering processes have 
been calculated only at threshold. They can be used to make universal predic- 
tions for cold atoms at sufficiently low temperatures. To determine the Efimov 
parameters and accurately and to test the accuracy of the universal pre- 
dictions, it is important to have universal predictions for rate constants as 
a function of temperature. This requires the calculation of the rates for the 
basic scattering processes as functions of the collision energy. The only such 
calculation thus far is the S-ware atom-dimer phase shift at collision energies 
up to the dimer break-up threshold [36] . For 3-body systems other than iden- 
tical bosons that also exhibit the Efimov effect, the universal results for the 
basic scattering processes have not even been calculated at threshold. 

The calculation of range corrections is important for improving the accuracy 
of universal predictions and also for quantifying their domain of validity. Al- 
though methods for calculating range corrections have been developed, they 
have not been implemented in a way that allows them to be easily applied 
to experiment. The calculation of 3-body observables in microscopic models 
for atoms near Feshbach resonance would also be useful for determining the 
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range of validity of the universal predictions. 

If the theoretical study of Efimov physics is still in its infancy, the experimental 
study of Efimov physics is a newborn. The beautiful results on ^^^Cs atoms 
from the Innsbruck group have provided the first glimpse of Efimov physics 
in cold atoms. These experiments have demonstrated that Efimov physics can 
have dramatic effects in alkali atoms in spite of their many deep 2-body bound 
states. They inspire confidence that Efimov physics will also be observable in 
other cold atom systems. 

The next milestone in the experimental study of Efimov physics in cold atoms 
will be quantitative tests of the correlations between 3-body observables pre- 
dicted by universality. For example, if a is large and positive, universality 
gives correlations between atom-dimer scattering and 3-body recombination. 
If a Feshbach resonance is used to tune the scattering length through ±00, 
universality gives correlations between 3-body recombination at large positive 
a on one side of the resonance and 3-body recombination at large negative a 
on the other side of the resonance. 

The ultimate confirmation of Efimov physics would be the experimental veri- 
fication of the discrete scaling symmetry. This could be accomplished by the 
observation of loss features in 3-body recombination or in atom-dimer scat- 
tering at two or more values of a that differ by powers of the discrete scaling 
factor. For 3-body recombination with identical bosons, the large size of the 
discrete scahng factor 22.7 and the scahng of the loss rate makes this diffi- 
cult. Atom-dimer scattering is more favorable in this regard because the loss 
rate scales only as a. Verification of the discrete scaling symmetry would be 
easier in a system consisting of two heavy atoms and a third light atom be- 
cause the discrete scaling factor can be much smaller than that for identical 
bosons. One particularly favorable possibility is the system consisting of two 
^'^'^Cs atoms and a ^Li atom near a Feshbach resonance in the Cs-Li interaction. 
This system has a discrete scahng factor of only 4.88 [86]. 

The discovery of the Efimov effect in 1970 provided a hint of the remarkable 
phenomena associated with Efimov physics that was waiting to be uncovered. 
For 35 years the study of Efimov physics was an almost exclusively theoretical 
enterprise that could be characterized as an exploratory mine shaft. The ad- 
vent of cold atom physics made it inevitable that Efimov physics would also be 
studied experimentally. This field offers a long list of systems in which Efimov 
physics could occur, as well as offering the possibihty of tuning the scattering 
length experimentally. An experimental mine shaft into Efimov physics has 
finally been opened by the Innsbruck experiment on ^'^^Cs atoms. This should 
stimulate the opening of mine shafts into Efimov physics in other cold atom 
systems. Extensive experimental effort to extract the ore and additional the- 
oretical effort to refine the ore will be required to reveal the full beauty of 
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Efimov physics. 
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